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1 Introduction 



The end of the Lep period represented a moment of transition in the development of techniques designed 
for producing high precision results for collider physics. 

It is a well known fact that the advent of a new hadronic machine, the LHC in this case, gives a privileged 
role to QCD but what is the real implication of this fact? 

QCD is a theory devoted to studying strong interactions by means of perturbative methods, a road 
which has been made possible by the discovery of asymptotic freedom. It is a theory with very few scales 
and therefore its perturbative aspects are technically simpler. For this reason most of the new ideas are, 
first of all, tested in QCD where we have an extensive literature of explicit results up to four-loop Feynman 
diagrams. The very recent twistor spinoff for collider physics [1] has its immediate target in the realm of 
QCD going beyond the traditional field-theoretic point of view. 

Generally speaking we are now witnessing the development of three parallel roads: techniques (mostly 
Monte Carlo driven [2] or twistor inspired [3]) to deal with tree level processes with many particles in the 
final state; complete one-loop calculations for 2 — > 4 processes with a particular emphasis on the proper 
treatment of intermediate unstable particles [4] and genuine electroweak two-loop calculations of physical 
observables ([5] and [6]), including supersymmetric effects [7]. The first item in the list is deeply linked to the 
LHC physic programme and has already collected a sizable number of important results. Extension of one- 
loop calculations to four particles in the final state represents a bridge between LHC and LC physics while 
two-loop electroweak physics is, at the moment, in some early stage of development. It is worth mentioning, 
however, that past history has told us about the importance of a complete calculation when dealing with 
the final analysis of the experimental data. 

Lepl physics has ben mostly dealing with 2^2 processes, where one has been able to assemble the most 
complete set of predictions in the whole history of radiative corrections. Most notable is the fact that the 
technology needed for the operation has been pushed well beyond what was expected at the beginning of 
the period, very much as in the case of the experimental analysis. The key ingredients have been 

- the construction of a complete one- loop renormalization procedure for the standard model [8], with 
the inclusion of several leading and next-to-leading higher order effects; 

- the development of libraries for assembling the relevant one-loop diagrams, including reduction of 
tensor integrals to scalar ones and analytical evaluation of the latter [9] ; 

- inclusion of higher-order, real and virtual, QED effects [10]; 

- development of fitting procedures to deal with pseudo and realistic experimental observables. 

At the very end of the Lep period it became evident that these procedures should be generalized if we 
want to have a full two-loop interpretation of the data. Furthermore, it became obvious that a step forward 
is necessary for treating the technological elements which are mandatory in obtaining a satisfactory solution 
to all the items of our wish-list. There are several aspects which deserve a specific comment: the main one 
can be summarized by saying that assembling a two-loop package is not, anymore, a one-man-show as it was 
in the past but requires instead a dedicated involvement of a well coordinated group. 

There are many new aspects in the project that put an unprecedented challenge: for any multi-scale 
process we have to abandon the fully analytical way and new, numerical, algorithms must be developed to 
handle the complexity of Feynman diagram evaluation with cutting-edge calculations. Generally speaking 
we are referring here to the evaluation of master integrals, a minimal set of diagrams to which all other 
diagrams are reduced. The techniques for reduction should be general enough to handle three or more 
external legs in a fully massive world. On top of that one should not forget to develop a comprehensive 
two-loop renormalization procedure, flexible enough to be implemented in practice. 

What we are describing here is a multi-step program and only one of the steps will be presented in 
this paper; to fully understand the results, however, it is important to mention that everything has been 
developed in parallel and what we are going to illustrate is already fitting in the general layout of our 
program. The general strategy for handling multi-loop, multi-leg Feynman diagrams was designed in [11] 
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and the collection of results necessary for evaluating two-loop, two-point integrals can be found in ref. [12]. 
Next the calculation of two-loop three-point scalar integrals: infrared-convergent configurations are discussed 
in [13] and infrared- and collinear-divergent ones are analyzed in [14]. Finally, our method for dealing with 
two-loop tensor integrals can be found in [15] and results for one-loop multi-leg integrals are shown in [16]. 

First of all, we needed a technique for reducing diagrams to some set of master integrals. Here, we are 
referring to a class of two-loop calculation with few external legs, avoiding the exponential complexity that 
would be encountered with many legs. This step requires our ability in treating all relevant limits, when 
masses are negligible or when external momenta must be set to zero. We should mention that the real world 
is massive, with all masses different, but most of them negligibly small. To mention an important aspect we 
may say that any numerical treatment that docs not extract collincar and Sudakov logarithms from the very 
beginning is doomed to failure. Furthermore, all the results must be stored in the most appropriate form, 
since our final numerical integration must be absolutely stable. 

Next, a full two-loop renormalization had to be designed and we have to make clear the meaning of 
this sentence. Diagrams have to be generated with some automatic procedure; in our project, we did not 
want to rely on any black box, so we have constructed our own set of procedures, creating the QraphShot 
package [17]. 

After assembling diagrams we wanted to perform all sort of canonical tests which, essentially, amounts 
to check all possible unrenormalized WST identities. After checking that everything has been properly 
assembled we have been moving to the next logical step, removing all ultraviolet infinities. At the very 
end one has to admit that the removal is a simple business but we pretend having done it in the most 
rigorous way. To mention an example we quote the problem of overlapping divergencies; non-local residues 
of ultraviolet poles must cancel in the total if unitarity of the theory has to be preserved: it is highly non 
trivial but they indeed do cancel. After that one has to check renormalized WST identities. 

Removing ultraviolet poles means trading bare parameters for renormalized ones and, after that, renor- 
malization equations must be written and solved which is equivalent to express theoretical predictions in 
terms of an input parameter set. In this paper, we illustrate the role of the running e.m. coupling constant 
and give a full account of the calculation of the Fermi coupling constant. 

Another, non trivial, aspect of our work is related to proving that what we expect to vanish is indeed 
zero; for instance one can prove that the standard model can be made, to some extent, QED-like: at the level 
of S-matrix elements we expect that vertices, with the inclusion of wave-function factors, do not contribute 
to the renormalization of the electric charge. To prove this property at two loops is far from trivial. 

It is by now clear that one cannot produce all the results in a two-loop calculation without having a 
software package, QraphShot [17] in our case; QraphShot collects all the relevant algorithms. At this stage 
the package is far from being user friendly (for instance we even miss a user's guide) and it is not even clear 
if somebody else will use it but writing the code has been an essential step without which none of the results 
presented in this paper would have been achieved. Once more, past experience is telling us that with the 
present level of complexity there will be no time, in a few years, to go back and to allow for extensions which 
were not foreseen from the very beginning. 

To summarize, we have performed several steps towards complete two-loop renormalization in the stan- 
dard model, steps that could be easily generalized, with few minor adjustments, to an arbitrary renormal- 
izable quantum field theory The introductory elements have been given in [18] (hereafter I) and in [19] 
(hereafter II). 

As it is well known, finite S'-matrix elements can be obtained without the explicit use of counterterms, a 
fact that has been fully described in the one-loop renormalization [10]. However, most of the people arc more 
familiar with the language of counterterms and, for this reason, we have decided to adapt our approach. 
In part II we have described the strategy for making the one- and two-loop Green functions of the theory 
ultraviolet finite; this amounts to rewrite the Lagrangian in terms of renormalized quantities. 

The next step in our procedure will be to express renormalized parameters in terms of physical observables 
belonging to some input parameter set (hereafter, IPS). In any one- loop calculation one tries to improve 
upon the accuracy of the result by including leading (and even next-to-leading) higher order effects (which 
are often available in analytical form) and by performing resummations which, quite often, cannot be fully 
justified. Here we focus on the issue of setting up a complete two-loop calculus, postponing the question of 
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building improved two-loop resummations. 

At two loops we have a new feature, as described in details by many authors: the use of on-shell mass 
renormalization is not allowed anymore and complex poles arc the only meaningful quantities: a property 
of the S -matrix, as it is often stated [20]. A better argument is that complex poles are gauge parameter 
independent to all orders, as shown by using Nielsen identities [21]. 

Renormalization with complex poles should not be confused with a simple recipe for the replacement 
of running widths with constant widths; there are many more ingredients in the scheme. Actually, this 
scheme allows for the introduction of a beautiful language, the one of effective (complex) couplings. The 
whole organization of loop corrections is most conveniently organized in terms of running couplings, as 
illustrated by the development of the fermion loop approximation. Unfortunately gauge invariance - the 
ingredient at the basis of the successful proof of renormalization - prevents us from fully extending the use 
of running, resummed, couplings to the bosonic sector of the theory (as a matter of fact even the concept of 
a fermionic sector is meaningless abeyond one loop). Despite this caveat we will organize the presentation 
of our results according to the language of running couplings; in a way this language is much easier and 
one should only remember that final results must be expanded in perturbation theory up to second order. 
Partial resummations, i.e. the attempt of isolating gauge parameter independent blocks that can be moved 
freely from numerators to denominators is beyond the scope of this paper. 

Neglecting the fine points of the procedure we may say that the transition from renormalized parameters 
to physical quantities depends on the choice of the latter. Quantities like G F and a, the Fermi coupling 
constant and the fine structure constant, will always be included in our choice of the IPS. For most of the cases 
under consideration fermion masses and the Higgs boson mass are only needed at one loop; here on-mass-shell 
masses can be used. For gauge bosons instead, we have to extract from the data some information about 
the position of the corresponding complex poles. The fact that we still indulge in presenting theoretical 
prediction for, say M W1 is a consequence of an established attitude that every pseudo-observable can be 
defined and derived from data, although its theoretical degree of purity is vanishing small. To an even larger 
extent the situation with the Higgs boson mass (as used in any LHC Monte Carlo program) is fully unclear. 

Having adopted this strategy, we may say that the final step in renormalization can still be seen as 
the moment where we write a system of (essentially) three coupled equations, relating g, sg and M (all 
renormalized quantities) to G F , a, .... 

In this paper, we shall follow the same notations as defined in our companion papers I and II. We therefore 
refer the reader to that papers for notations; in particular, M stands for the bare (renormalized) W boson 
mass (we do not distinguish unless strictly needed), Mq = M/cg where eg is the bare (renormalized) cosine 
of the weak-mixing angle; M w is the on-shell W boson mass and M z = M w /cg. 

The outline of the paper is as follows: in Section 2 we summarize our procedure, and in Section 3 
we present a renormalization equation based on the use of the Fermi coupling constant; aspects of the 
calculation connected with the proper definition of 7 5 are discussed in Section 4. A second renormalization 
equation connected with the fine structure constant is analyzed in Section 5. The running of a beyond 
one-loop fermion terms is discussed in Section 6. Complex poles are introduced in Section 7 and solutions 
of the renormalization equations in Section 8. Loop diagrams with dressed propagators are introduced in 
Section 9, unitarity, gauge parameter independence and WST identities in Section 10. A suggestion on how 
to improve the complex mass scheme [4] is finally introduced in Section 11. 

2 Outline of the calculation 

The whole renormalization procedure has been summarized in the flowchart of Fig. 1 where IPS stands 
for Input Parameter Set. Object of this paper is to introduce renormalization equations (RE) and to solve 
them, therefore undertaking the task of writing any (pseudo-)observable, not in the IPS, in terms of the 
quantities of the IPS. The minimal standard model is essentially a three-parameter theory and therefore we 
seek for a system of three (coupled) REs whose unknowns are g,M and sg, all renormalized parameters. 
In Section 3 we discuss the first RE, related to the Fermi coupling constant, G F . In Section 5 we present 
the second RE, related to fine structure constant a. In Section 7 we introduce the third RE, based on the 
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Figure 1: Renormalization - flowchart. Feynman rules define the theory, renormalizability guarantees that ultraviolet 
poles have polynomial residues and can be subtracted. Any input parameter set allows us to replace renormalized 
quantities with experimental data. A prediction follows. 



notion of complex pole for unstable gauge bosons. The solution of our RE is discussed in Section 8 where 
we present two choices for the IPS. 



3 The Fermi-coupling constant 

In this section we present the results for G F . A critical observation is that writing a renormalization 
equation for the Fermi-coupling constant should not be confused with predicting the muon lifetime or the 
Fermi coupling constant itself (also in an effective field-theory approach, see i.e. Ref. [22]). 

To proceed further, we illustrate our method for relating the A/S'-renormalized parameters of the Standard 
Model to the Fermi coupling constant at two loops and we construct our first renormalization equation. 

The Fermi coupling constant, G Fl is defined by 

1 Gl ml 

where t m and m^ are the observed muon lifetime and mass. The parameter Aq summarizes both real and 
virtual corrections to r M , at leading order in G F and to all orders in the fine-structure constant a. The 
corrections are generated by the effective Lagrangian 

G 

Ajedxft = £ft + £qed, £ ft = ~y= [tV7a (1 + 75) M] ' [eja (1 + 75 ) V e ] . (2) 

Here C FT describes the contact-interaction Fermi theory (FT), £ QBD is the usual QED Lagrangian and fi, e, 
and v e are the spinor fields for the muon, the electron and their related neutrinos. 

The inclusion of electromagnetic effects allows to determine the numerical value of G F in terms of the 
measurable quantities r M , a, m M and the electron mass m e . The combination of the one- loop [23] and two- loop 
[24] QED x FT contributions with the non-perturbative hadronic components [25] leads to estimate 

G F = (1.16637 ± 0.00001) x lO^GcV" 2 . (3) 
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Unfortunately, another definition is often employed in the literature (see e.g. Ref. [26]), 



where M w is the W-boson mass and F{m 2 e /m 2 l ) follows from the phase-space integration at lowest order in 
perturbation theory, 

F{x) = 1-%x-12x 2 In a; + 8 a; 3 -a; 4 . (5) 

However, the factor (1 + 3m 2 ( )/(5M^) is the tree-level VF-propagator effect and is not generated by the 
Fermi-contact interaction. Moreover, as pointed out by the authors of Ref. [24], the function F does not 
factorize in the same way at higher orders. Therefore, in order to avoid unnecessary ambiguities, we will use 
the definition of Eq.(l). 

In the context of our renormalization procedure, G F is an input data and we have to derive the appropriate 
renormalization equation. First we express t m through the SM renormalized parameters. Next, we match 
our result with the definition of Eq.(l) getting a relation between a measurable quantity, G F , and the SM 
renormalized parameters, 

Here g and M are the renormalized weak-coupling constant and TU-boson mass and the quantity Ag is 
constructed order-by-order through the purely-weak corrections to the muon lifetime. Electromagnetic com- 
ponents are already included in Eq.(l) and they are discarded through the matching procedure. 

It is worth noting that the answer for Ag will contain renormalized parameters and counterterms. The 
choice of the renormalization scheme determines the explicit expressions for the counterterms and the final 
result for Ag. A popular strategy, followed in recent two-loop calculations of the muon-decay width [27,28], 
is to employ the on-mass-shell (OMS) scheme and to define renormalized parameters by means of measurable 
quantities. At lowest order one has simple relations, 

M 2 

e 2 =g 2 s 2 e =47ra 1 c 2 = M = M W . (7) 

z 

Here e is the renormalized electric charge, M w and M z are the on-shell masses of the W and Z bosons and 
cg(sg) is the renormalized cosine (sine) of the weak-mixing angle. A replacement of the OMS renormalization 
prescription into Eq.(6) leads to the traditional parametrization employed to describe the interdependence 
between the masses of the vector bosons introduced in Ref. [29] . 

Rather that defining exactly what the renormalized parameters are, like in Eq.(7), we prefer to follow 
a minimal-subtraction scheme. Let us prescribe the values of the counterterms as an intermediate step to 
remove ultraviolet poles and regularization-dependent factors. Relations among renormalized parameters and 
physical quantities are not imposed by hand and represent the solution of the chosen set of renormalization 
equations. Therefore, our result for Ag should not be confused with Ar as reported by the authors of 
Ref. [27,28]. 

The basic prerequisite to derive the renormalization equation for G F is the extraction of the QED x FT 
components from the full SM calculation of the muon lifetime. In Subsection 3.1 we discuss our method and 
detail the result with some relevant diagrammatic examples. Note that the definition of G F relegates the 
infrared structure of muon decay to the soft electromagnetic effects summarized by Aq in Eq.(l) and does 
not affect the hard weak remainder, Ag. 

Moreover, in our approach, external legs will be provided with appropriate wave-function renormalization 
(WFR) factors; since most of the existing literature deals with the OMS scheme, where WFR factors are 
usually replaced by field counterterms, we devote Subsection 3.2 to discuss their role. Of course, the connec- 
tion between field renormalization and wave-function factors is well understood, they are simply connected 
by a field transformation. 

Finally, in Subsection 3.3 and Subsection 3.4, we give our results for Ag at one loop (consistency check) 
and at two loops. 
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Figure 2: Diagrammatic interpretation of Eq.(13). The first graph is a box diagram in the full SM context. The 
one-photon vertex diagram is a QED correction in the Fermi-contact interaction, denoted by the black circle. 



3.1 Extraction of the electromagnetic components 

In this subsection we consider the Standard Model tree-level amplitude for muon-decay and fix our 
notations an conventions. The process is 

-> iV (Pa) + e(p 3 ) + V e (p A ). (8) 

Neglecting the electron mass and the squared momentum carried by the intermediate VF-boson propagator, 
the result in 't Hooft-Feynman gauge reads as 

M° = - (2ir) 4 i ^ [u 2 7q (1 + 7s) ui ] ■ [u 3 la (1 + 75) Vi] , (9) 

where we used the short- hand notation for spinors Ui — u(pi) and Vi = v(pi). 

The key observation is that the QED corrections in the context of the Fermi-contact interaction can be 
systematically identified and removed from the full SM amplitude. The one-loop matrix element factorizes 
as 

M 1L = M° ( 1 + AAf ^ DXFT ) ( 1 + AMl L EAK ) + O {g 6 ) , (10) 

where AA4q EDXFT is a scalar function which summarizes all the QED x FT soft contributions and the purely- 
weak remainder is entirely relegated to the hard component A.M^, AK . At two loops the amplitude can be 
decomposed as 

M 2L = M° ( 1 + AMl L EDXFT + AM 2 Q L EDXFT ) ( 1 + AMl L EAK + AM 2 W L EAK ) + O ( 5 8 ) , (11) 

where AA^qedxft an d ^A^widak denote two-loop soft and hard corrections. The class of soft effects generates 
the parameter Aq of Eq.(l). We discard them and we identify the hard terms with the quantity Ag introduced 
in Eq.(6), 

Ag = AMl L EAK + AM 2 W L EAK + O (g 6 ) . (12) 



A simple one-loop example shows how the QED x FT components factorize in the SM amplitude. We 
apply the procedure introduced in Ref. [30] to separate off the soft electromagnetic corrections from the hard 
remainder in the SM box diagram of Fig. 2. If q denotes the momentum of the virtual photon and q +p the 
momentum flowing through the M^-boson line, we decompose the integrand in the large PF-mass limit, 



q 2 [{q+ p) 2 + M 2 } M 2 \q 2 q 2 + M 2 J \ q 2 + M 2 

Here the first term generates the soft QED-vertex correction to the local Fermi interaction of Fig. 2, which 
can be obtained replacing the T^-boson propagator by 1/M 2 . We discard this component and we evaluate 
the hard weak remainder in the large vector-boson mass limit, where the external momenta and the lepton 
masses can be safely neglected, 

4 2 r 

WEAK 8M2 ^ J a 9( g 2)2( 9 2 +M 2) i ^- ^) 



G 



Figure 3: Two- loop box diagrams for muon decay. 



Here [i is the 't Hooft unit of mass and the spinor chain reads as 

= [u2 7a7+7M 7/3 U l] • [«3 7/3 7^7q7+W4] , (15) 

where j± = 1 ± 7 s . Note that the purely- weak components can be obviously obtained starting from the full 
SM amplitude and subtracting the QED x FT terms. Therefore, instead of extracting the hard correction 
from the box diagram of Fig. 2, we can use an alternative strategy: 

I) We construct the difference of the box diagram and the QED x FT component of Fig. 2. This can be 
evaluated nullifying the external momenta and the lepton masses. 

II) Here the soft electromagnetic term is a massless tadpole which vanishes in dimensional regularization. 
With a non-zero electron mass, of course, it develops infrared and collinear singularities. Therefore, any 
information about the infrared structure is lost but these effects are already included in Aq of Eq.(l) and 
are not relevant for Ag. 

III) As a result, the hard part follows directly from the complete one-loop diagram by nullifying the lepton 
masses and the external momenta, 

4 2 p 

M wj = _9 S B A-n d n pW 7 

WEAK g A* J a 1 ( g 2)3( g 2 + M 2) ■ 

Using the fact that in dimensional regularization an integral without scales is zero it is easy to see that the 
two representations of Eq.(14) and Eq.(16) are completely equivalent. Since QED factorization should not be 
proven by assumption, and one has anyway to separate off the QED x FT components, the choice between 
the two procedures is just a matter of taste. However, the subtraction method of Eq.(16), introduced in 
Ref. [28], appears more appropriate in the context of an automatic approach, where all the Feynman diagrams 
are generated and evaluated by neglecting the soft scales from the very beginning. 

The matching on the Fermi-theory spinor-chain configuration can be finally completed when we introduce 
the projector, 

V = 1^4 7a 7+ "2] • [ui7a7+m 3 ] (17) 

spins 

which will act on the matrix element for muon decay. For example, for the hard part of Eq.(14) we can write 
M W1 =M° ■ AM W1 with AM "~ 7 = WEAK (18) 

WEAK J 1 J 'WEAK! vviuii "J"' WEAK 'PMl ' 

Two-loop virtual corrections to the muon-decay amplitude can be classified according to Fig. 3. Here 
diagram a) contains two photons and one hard scale and can be treated in complete analogy with the one- 
loop case. After shrinking the heavy line to a point like in Fig. 4, the two-photon graph in the local Fermi 
theory is discarded. The hard weak remainder is obtained by nullifying the soft scales in the full two-loop 
box diagram. Diagram b) includes just heavy components, and does not require any subtraction. We will 
evaluate it in the soft limit obtaining a finite answer. 

Since in diagram c) soft and hard components are entangled, we discuss it in more detail. We consider 
the subloop with a photon and a Z boson and we decompose it in the large Z-boson mass limit. If q\ is the 
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Figure 4: Two-loop soft ® soft splitting. 
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Figure 5: Two-loop soft <g> hard splitting. 



momentum flowing along the photon line and q\ — (72 is the momentum of the Z propagator, we can write 



2 



1 1 1 ( . , 2g x • q 2 - q( 



1+ 2 ;%^ +... , (19) 



9? [(ft-? 2 ) 2 +M 2] g 2 <z 2 + M 2 V q$ + M 2 

where Mq — M/cg. The first term gives the product of a one- loop QED x FT vertex and a one- loop weak 
remainder, as shown in Fig. 5. The rest is a purely-weak two-loop component which will be evaluated after 
nullification of external momenta and lepton masses. 

The soft x hard term is crucial in proving the two-loop factorization property of Eq.(ll), because it 
generates the AM\^. UXWT x AM^ EAK part. To evaluate the hard remainder, we start again from the 
complete representation for the SM diagram, where we nullify the soft scales, 

kytWZ-y _ _ 1 J_ 9 2 s l ,.2(4-n) f m m QlS l2v glp Q2a pWZ 7 /r,^ 

weak - (2?r)4z 32 Cg J qid q2 (g 2 ) 3 fe 2 ) 2 fe 2 + M 2 )[(gi-«z 2 ) 2 + M2] L '"p°- [Z[)) 

Here the spinor chain reads as 

T ™vpa = [U2 7« 7+ 7^ 7/3 (v» + 7« 7t u l ] ' [ "3 7>y 7p 7/3 («e + OeTs) 7a 7« 7+ ^4 ] (21) 

and u/ = /j 3 ** — 2Q 2 Sg and a/ = where /j 3 ** is the weak-isospin third component for the fermion / 
and Qf is the related electric-charge quantum number. Finally, a correct description of muon decay involves 
also real radiative effects. This is essential in showing that the QED x FT corrections are finite [23,24] and 
allows to define the Fermi-coupling constant through Eq.(l). Nevertheless, one can easily prove that real 
soft-radiation diagrams can be systematically neglected from the calculation. 

3.2 Wave-function renormalization factors 

Wave-function renormalization factors for external legs (hereafter WFR) contain derivatives of two-point 
Green's functions which are infrared divergent and deserve a special discussion. In order to derive explicit 
representations for the muon-decay case we consider the fermion one-particle irreducible Green's function, 

E(p) = (2ir) 4 i (a s (p 2 ) + [a v (p 2 ) ~ a A (p 2 ) l5 ] if}, (22) 

where we introduced a scalar, a vector and an axial form factor (the SM pseudo-scalar component obviously 
vanishes). The fermion Dyson-resummed propagator reads as 

S(p) = {(2n)U [(l - a v {p 2 ) + a A {p 2 ) lb yi> + m - a s (p 2 )] }~\ (23) 
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where m is the renormalized fermion mass. Let us accordingly define appropriate WFR factors using a 
second representation for S(p) around the on-shell fermion mass, m os , 

S(p) = {(2tt) 4 i [( 1 - W v + W A 75 ) (if + m os + O {{if + m os ) 2 )) ( 1 - W v - W A l5 )] }~\ (24) 

Here W v and W A are the vector and axial WFR factors, whose expressions can be obtained through fermion- 
mass renormalization followed by a straightforward matching procedure of Eq.(24) with Eq.(23). 

Fermion-mass renormalization . We trade the renormalized mass for the on-shell one writing 

1 



S( P )D(p z ) = - 



{2ir) 4 i 



1 - a v (p ) + a A (p )7 5 if - m + a s (p ), 



(25) 



D(p 2 ) = [1 - 2a v ( P 2 ) + a 2 v (p 2 ) - a 2 A (p 2 )} p 2 + [m - a s (p 2 )} 



and imposing fermion-mass renormalization, D(—m 2 os ) — 0. We obtain a pcrturbativc solution expanding 
a s , a v and a A through powers of the weak-coupling constant g, 



aAp 2 ) - ]T.9 2 XV), I = S,V,A, 
i=i 

and performing a second expansion around the mass shell, 

a] L {p 2 ) = A\ L + 2m os {if + m os ) B\ L + O ({if 



+ m c 



) 2 )- 



(26) 



(27) 



Here we introduced short-hand notations for the form factors evaluated on the fermion mass shell and their 
derivatives respect to the squared external momentum, 



A] L = a] L {-m 2 os ), B\ L 

The solution for fermion-mass renormalization reads as 

9 2 (A\ L -m os A v L ) +9*{A¥ 



da] L {p 2 ) 



dp 2 



m 



A 2 V L + \(A\ L ) 2 ]}+0{ 9 «) 



(28) 



(29) 



and removes from Eq.(23) the renormalized mass. Note that one- loop fermion-mass renormalization involves 
just the axial and scalar form factors evaluated on the fermion mass shell. At two loops, instead, we are left 
to consider an additional term given by the square of the one-loop axial component. 

Wave-function renormalization. Expand the WFR factors introduced in Eq.(24) order- by-order in pertur- 
bation theory, 

oo 

Wj = 9 2j W i L ' I = V,A. (30) 

We wish to determine the explicit expressions of the WFR factors at one and two loops by matching Eq.(23) 
with Eq.(24), 



VF 



1L 



IL 

os u s 



m os B v , 



W] L = ^A 



1 



1L 



VF, 



2L 



= - I A 2 V L 
2 I v 



(A?) + (A\ L ) ]} + m os (Bf + tA v L B\ L 



2 

m os 



W 



2L 



{Bf + \ [Ai L Bi L (Bl L ) 2 ] }- m* oa B?B? + \m% s ( Bv L )\ 



Af + -Ai L A\ L 



- m„ s A\ L Bl L 



— m„ s A, B, r . 



(31) 
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Beyond the tree level we provide the external spinors with the WFR factors. The one-loop matrix element 
includes a component given by tree level diagrams where each external leg gets separately a WFR factor at 

g 2 (W^ L + W A L l5 ) u, g 2 u ( W X V L - W A L l5 ) . (32) 

The two-loop amplitude requires a more careful treatment, and involves tree classes of diagrams containing 
WFR factors: 

- tree level diagrams where three external legs are not corrected, and one includes a WFR factor at 
0(9'), 



•s 4 



4 — 

■9 u 



75 



W 2L + ( W < L ) 2 + { W a L + 2 Wl L W\ 

I = V,A 

w 2L + J2 ( W ' L ) 2 - ( w a l + ^wl L w\ L ) l5 



u, 

(33) 



I = V ,A 



- Tree level diagrams where two external legs are simultaneously providedwith a WFR factor at O (g 2 ) 
as in Eq.(32). 

- One-loop diagrams where the external legs are corrected one-by-one at O (g 2 ) through the WFR 
components of Eq.(32). 

Concerning the muon-decay case, we are dealing with an external positron. Relations for antiparticlc 
spinors follow trivially from Eqs. (32)-(33) replacing u and u by v and v. 

Eq.(31) contains derivatives of the 1PI Green's function. These terms are suppressed by positive powers 
of the fermion mass, and in QED they survive with double and single infrared poles and logarithmic mass 
singularities. In the computation of the purely-weak component, instead, we can neglect the fermion masses 
in Eq.(31) and compute WFR factors through the scalar, vector and axial form factors evaluated at zero- 
momentum transfer. 

3.3 Two-Loop Corrections to Ag 

In this section we will extend to the two-loop amplitude; in order to evaluate Ag we follow three steps: 

- We generate the one-loop and two-loop matrix elements A4 1L and Ai 2L for muon decay by nullifying 
the external momenta and the light-fermion masses. As discussed in Section 3.1, this implies an 
automatic subtraction of the soft QED x FT corrections. Therefore, our expressions for the one-loop 
and two-loop amplitude only contain the weak contributions. 

- We project A4 1L and Ai 2L on the tree level amplitude through the projector defined in Eq.(17) and we 
extract the hard components AM^ EAK and A.M^, AK . The sum over the spins of the external particles 
leads to traces over Dirac matrices which have to be evaluated in n dimensions. 

- At this stage, we have to compute a large set of massive tadpole diagrams. Using integration- by-part 
identities they are reduced to one master tadpole integral, as shown in appendix A of II. 

Before showing the two-loop result, we review the one-loop corrections. The one-loop hard contributions 
can then be conveniently written as 

AM 11 — AM 1L +AA// 1L f34l 

OJ,l WEAK " J l WEAK,W ' WEAK, REST 5 V°^V 

where AMJ^ak.w denotes universal H^-boson self-energy corrections, and AA4^ EAKtREST represents process- 
dependent vertex, box and wave-function components. The result for the latter reads as 



2 



AM 1L = — 

WEAK. REST 1e 

lb 7t 



7-4 *2 
2s 2 



Incj . (35) 
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Since we removed all ultraviolet poles and regularization-dependent factors from one- loop diagrams, the 



quantities AM 



1L 

WEAK,W 



and AMhh 



WEAK. REST 



do not show any ultraviolet-divergent component. However, as 



explained in [31], AM^ EAK REST is ultraviolet finite also in the bare theory, owing to the introduction of the 



parameter Ti (see I). 

Two- loop corrections to the muon-decay amplitude can be organized through a universal component, 
represented by T^-boson self-energy reducible and irreducible diagrams, and a process-dependent part con- 
taining irreducible two-loop box, vertex and WFR diagrams and reducible ones, 



AM 2L 

'weak 



Irreducible two-loop box diagrams give (xi 



AM 2L 

m 2 /M 2 ) 



WEAK,W 

2 /A/f2\ 



AM 



AM 



2L 

WEAK, BOX 



UV 



9 



(16tt 2 ) : 



2L 

WEAK, REST " 



157 3 1 



(36) 



12 



+ 7— + 

4d 



157 



83 



- %4xt - — (-T - 12^ + 2 + 12x 2 t In x t + \ In c 2 ) 

Z X H \Cg Cg / 

-^x H s 2 - (6 - ~ - 4c 2 . + 2cfj lnc 2 



(fin*. 



lnxa + 3x t sl lnx t \ 

z H - 1 J - 1 

and the sum of the other process-dependent components is 



2 Se x t 



A M 2L 

weak, other 



= - AM 



2L 

WEAK, BOX 



(37) 



(38) 



where we have included one-loop diagrams with a one-loop counterterm insertion but no two-loop countert- 
erms, since one can show that - exactly like for the one-loop case - their contribution cancels, and the sum 
of Eqs. (37)-(38) is an ultraviolet-finite quantity by itself. Therefore, a one-loop renormalization is enough 
to cancel the ultraviolet poles for A_M^. AK REST and to construct the process-dependent component of the 
Fermi-coupling constant by neglecting two-loop counterterms. Note that a crucial role here is played by 
the parameter IV The inpact of two- loop corrections on the finite quantity 8 G is shown in Tab. 1. The 
complete expressions of all the components of Ag are too lengthy to be presented here and have been stored 
in http: / /www. to. infn.it / ~ giampier/REN / GF.log. 

3.4 Process independent, resummed, Fermi constant 

If we neglect, for the moment, issues related to gauge parameter independence it is convenient to define 
a G constant that is totally process independent, 



Ag = 6 G + Ag s 



G 



Gj_ 

V2 



1 - 



;m 2 



60 



E 



9 



16 71- 



x(") 

°G ■ 



(39) 



Alternatively, but always neglecting issues related to gauge parameter independence, we could resum 5 G by 
defining V2G R = G F /(1 + S G ). In one case we obtain 



G = 



8M 2 



1 - 



16tt 2 M 2 



,(0) 



J ww(C) — £^^(0) + 



,(2) 



16 7T 2 



(0), 



where Y. ww is the W self-energy, whereas with resummation we get 



G R = 



9 



8M 2 



1 



9 



16 7T 2 M 2 



,(0) 



,(0) 



.(i) 



(0) 



16 7T 2 



^(2) 



(0) 



^(2) 



(0)<5. 



(i) 



(40) 



(41) 



In Section 8 we will show how solutions for the renormalization equation of the SM are simpler when written 
in terms of G. 
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4 The 7 5 problem 



To compute two-loop corrections to fi decay we must specify how to treat 7 s . The naive scheme is defined 

by 

{7^ , 7 5 } = 0, Tr ( 7 5 7^ Y 7 Q 7") = 4 e^ af3 , (42) 

in 11— dimensions. A complete set of calculational rules with algebraic consistency requires to consider 7^, 7 
and e^ ua f3 as formal objects[32] with the following properties: 

{7", 7"} = 2$"" J, Trl = 4, 

7 5 = \ { ^ Tr (7^7*7/3) = ^ e^" 3 Tr ( 7 £7»7^) , 
{7^, 7 5 } = 0, [f, 7 5 ] = => A» = {7", 7 5 } = 2 7V, 
Tr 7 5 = Tr ( 7 5 7|t7l/ ) = Tr ( 7 5 7/n . . . 7/12p+1 ) = 0, 

Tr (7 5 7'V7 a 7 /3 ) = 4e^ Q/3 ( 7 5 ) 2 = 1. (43) 

Eq.(43) defines the 't Hooft - Veltman - Brcitcnlohncr - Mason scheme (hereafter HVBM). 

The HVBM scheme breaks all WST identities (so-called spurious or avoidable violations) which can 
be restored afterwards by introducing suitable ultraviolet finite countcrtcrms. The procedure, however is 
lengthy and cumbersome. 

The usual statement that we find in the literature is: consider the two diagrams of Fig. 6, they are 
the only place where the difference between the two schemes is relevant. If we can prove that the fermion 
triangles inside the diagrams of Fig. 6 in the HVBM scheme give 



1 ssW 



C ' t - — 1 ssW 



+ O (n - 4) , (44) 

ac 7 5 



(and for both schemes the anomaly is correctly reproduced) then the difference (0(n — 4)) is irrelevant 
since the remaining integration in the two-loop integrals is ultraviolet finite; the latter follows in general 
from renormalizability and it is confirmed by our explicit calculation which shows that all e— tensor terms 
are ultraviolet finite. However, Feynman rules are derived from the Lagrangian and one cannot change the 




Figure 6: Example of two-loop diagrams contributing to the anomaly. 

former but not the latter. The relevant analysis has been performed by F. Jegerlehner [33] and we simply 
repeat the argument; in constructing the SM Lagrangian we use chiral fields and the relation 

iPl,r 7 M ^r,l = 7j4>J^l-ip, (45) 
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is valid only if {7^ , 7 s } = 0. The consequence is that we miss a chirally invariant dimensional regularization, 
i.e. we cannot regularize fcrmion-loop integrals by continuation in n. If we insist on the trace condition then 
gauge invariance must be broken in order to obtain a pseudo-regularization, i.e. we use an n-dimensional 
fermion propagator 

° - - lf> + m P 2 =f+f. (46) 



F p 2 + m 2 — i5' 

Pseudo-regularization introduces spurious (i.e. avoidable) violations of WST identities which must be re- 
stored afterwards by introducing suitable ultraviolet finite counterterms. The procedure is lengthy because 
all Green functions containing a fermion loop are now anomalous and not only the fermion triangles inside 
the diagrams of Fig. 6. A simple example comes from the difference between one-loop naive and one-loop 
pseudo-regularized vector - vector, vector - scalar and scalar - scalar transitions; define finite counterterms 
in the pseudo-regularized formulation that make this difference zero. We will get a set of equations, 





V N c f £jf = 0. 



/ 

*r,ww = gf^ ^ [i (! + 2 W 2 + «f + m %)\ . ( 4? ) 

etc, where is the fermion color factor and where we have introduced the evanescent operator 

p 2 = (4-n)£ p f. (48) 



If we expand finite counterterms, 



etc, it is easy to derive a solution, 



\ + J—8Z w , (49) 
16 n z 



6Z W = ^(1 + 2E P ), SZ Z = ^ (^-8s 2 e -2^ (1 + 2^,), (50) 

etc. Although a solution for self-energies can be obtained it is clear that the number of anomalous WSTI is 
greater than the number of Lagrangian counterterms, i.e. counterterms specifically associated to parameters 
and fields. Even this fact does not pose a serious problem since we are talking about finite counterterms 
and, in principle, we can associate an ad hoc counterterm to each avoidable anomaly. Actually there is 
more; if all fermion loops induce anomalies then the argument of Eq.(44) is violated and we must introduce 
counterterms at O (n — 4) . 

Clearly it is not an ideal solution. After considerable wrangling one is lead to the conclusion that the 
only sensible solution is the one proposed by Jegerlehner [33]: n-dimcnsional 7-algebra with strictly anti- 
commuting 7 s together with 4-dimcnsional treatment of the hard anomalies. 



5 The fine structure constant 

To discuss the effect of radiative corrections on a and its interplay with renormalization we need the 
photon propagator, including two-loop contributions; using the results obtained in Sect. 5 of I, we write 



D AA (P 2 ) = 4 Y, 



n=l,2 



9 



16 7T 2 



n 



(n) 

QQ ; ext 



(P 2 )P 2 



(51) 



where the external LQ decomposition of self-energies in the neutral sector has been introduced and discussed 
in Sect. 6 of I (Eqs.(123)-(125)). Using e 2 = g 2 s 2 we write 



IWP 2 ) 



16 TT 2 



n^ ;oxt b 2 ,{™}) 



n 



(2) 



lg7r 2 s 2 QQ ; cxt 



(p 2 , {TO} , 4) 



(52) 
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where {m} denotes the full set of masses, including bosons; note that the two-loop SM contribution to the 
photon self-energy is not proportional to e 4 ; as a consequence, charge renormalization does not decouple 
from the remaining REs. Finally, we construct n QQ (0), as illustrated in II; for this operation several steps 
are needed: although one could start from the very beginning with p 2 = where all diagrams are vacuum 
bubbles that can be reduced by integration- by-part techniques [34] (but see also ref. [35]) to one master 
integral or products of one- loop functions, we prefer to perform the limit p 2 — > of the output of QraphShot 
which represents a highly non-trivial test of the procedure. It is useful to define 

ff(2) -_L n (2) (*>X\ 

11 QQ;ext — 2 QQ;ext' V JO ) 

and our renormalization equation reads as follows: 

^(0)^ + ^ + 7^:1^-1 = 0, (54) 



16 7r 2 Ana 



where a is the fine structure constant. Strictly speaking, our renormalization equations form a set of coupled 
equations; at one loop there is no residual dependence on the weak- mixing angle once we write e 2 = g 2 Sg, 
but the two-loop contribution modifies this simple structure. In this case, however, we should simply insert 
the lowest order result for s^-dependent terms of O (e 4 ) in Eq.(54). What to choose depends on the IPS 
that we select; postponing, for a moment, this question we obtain 



e 2 _ ^ a 



n ««(°) + ^ { [n«(o)] 2 - n^ ;cxt (o, s - 2 )}. (55) 



An a Air 

The value of s 2 in the argument of the two-loop corrections is fixed by the corresponding lowest order 
solution, for which we will need the remaining two renormalization equations. 
The total contribution to vacuum polarization is split into several components: 

n w (o) = n£(o) + n&(o) + ng-(o) + n^», (56) 

where the fermionic part contains three lepton generations, a perturbative quark contribution and a non- 
perturbative one. The latter, associated to diagrams where a light quark couple to a photon, is related to 
^ a had(-^z) [36]- The top contribution at one loop and two loops can be computed in perturbation theory 
due to the large scale of the top mass; always at two loops, diagrams where quarks are coupled internally to 
vector bosons are also computed perturbatively. However, QED and QCD contributions to the light-quark 
part are always subtracted. 

There are important tests on the result once the limit p 2 — > has been taken. QED is always included 
in the leptonic sector and collincar logarithms arc present in the final answer. However, as it is well-known, 
they are of sub- leading nature, i.e. the correction factor is proportional to g A In to 2 and not to <? 4 ln 2 m 2 ; 
indeed, the leading logarithms are controlled by the renormalization group equations and are related to 
Dyson re-summation of one-loop diagrams. 

For the quark sector we introduce a perturbative contribution and a non-perturbative one associated to 
diagrams where a light quark couple to a photon. Consider a doublet of light quarks, we subtract the QED 
part from the total and give the perturbative component; this means that for each up and down quark we 
neglect the two loop-diagrams which are built with all possible insertions of a photon line in a loop of light 
quarks and we also neglect the QED component in the light quark mass renormalization. It is worth noting 
that the subtracted terms are gauge invariant. Since QED has been subtracted we do not expect collinear 
logarithms. 

We also have a perturbative heavy-light contribution where the up quark is the top. In this case only 
the QED component of the 6-quark is subtracted and we do not neglect the top mass. Note that logarithms 
of the mass of the b quark are expected to disappear from the answer, as for any other light quark. In the 
following section we sketch how the result is constructed in QED. 
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5.1 The QED case 

To understand renormalization at the two-loop level we consider first the case of pure QED where we 
have 

2 4 

U QED (s, m) = ^ m) + n£2 



(57) 



where p 2 = — s and where we have indicated a dependence of the result on the (bare) electron mass. Suppose 
that we compute the two- loop contribution (3 diagrams) in the limit m — 0. The result is 



n^M) 



7 + 



(58) 



where n = 4 — e. This is a well-known result which shows the cancellation of the double ultraviolet pole as 
well as of any non-local residue. The latter result is related to the fact that the four one-loop diagrams with 
one-loop counterterms cancel due to a Ward identity. Let us repeat the calculation with a non-zero electron 
mass; after scalarization of the result we consider the ultraviolet divergent parts of the various diagrams. 
Collecting all the terms we obtain 



ITqq (s,m) 



4 1 + 24 ■ 



192 



1 



In 



/?(M 2 ) + 1 



s 2 /?(/z 2 ) /3( M 2 ) 



•0(1). 



(59) 



Note that the m dependent part is not only finite but also zero in the limit s — > 0; indeed, in the limit s — > 
and with /i 2 = m 2 / s — i S we have 



/3=(l-4 M 2 ) 



1/2 



j-5 + O (m„-3) 



so that 



n$(o,m) = -- + n 



(2) 
QQ ; fin 



(0,m). 



(60) 



(61) 



Eq.(61) is the main ingredient to build our renormalization equation and contains only bare parameters, 
respecting the true spirit of the renormalization equations that express a measurable input (a in this case) 
in terms of bare parameters (e and m in this case) and of ultraviolet singularities (this last aspect can be 
avoided by introducing counterterms). 

To make a prediction, the running of a in this case, is a different issue and, actually, does not depend on 
the introduction of counterterms: the scattering of two charged particles is proportional to 



l + f(s) + f(s) + 



W(*) 

Renormalization amounts to substituting 



16 7T 2 



(*) + 



n$( S ) + 0(e 6 ). (62) 



(16 7T 2 ) 2 



4 7T a -a 2 n«(0) 



|^{[nW(o)] 2 -n^(o)} + o( a 4 ), 



with the following result: 



= 4 * a { i + £- n« (s) + (£) 2 [n« W n« W + n?> W + o (« 3 ) } , 



!-/(*) 

nW(a) = nW(a)-nW(o) 



(63) 



(64) 



To have an ultraviolet finite result we also need that the poles in IIqq(s) should not depend on the scale 
s. This is obviously true for the one-loop result but what is the origin of the scale-dependent extra term in 
Eq.(59)? One should take into account that 



n^(s,m) = -~ + ^ ln^- + (1 + 2— )/?( — ) In 



1 



m 2 [3{m 2 1 ' s) + 1 



3 e 



M 2 



(3{m 2 /s) - 1 



20 4 , 

-y + 3 A ^ 
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and that m is the bare electron mass. To proceed further we introduce a renormalized electron mass (not to 
be confused with the physical one) which is given by 



m = m F 



16 7T 2 



6 



finite partj 



(66) 



If we write m 2 = m 2 {1 + 8) then 



s s 0{-*)8 

ln p(m 2 /s) + 1 = ^ f3(m 2 Js) + 1 



5 + 0(S 2 ) 
6 



P(m 2 /s)-l P(m 2 R /s)-l f3(m R ) 
Inserting this expansion into our results we obtain 



o(s 2 



(67) 



II QED (s,m R ) 



12 



5 1 m 
36 ~ 3 ~s~ 



2 , 



+ 



showing cancellation of the ultraviolet poles in II„ '(s, m R ) with n = 1, 2. 



1 , m 2 „ 1/1 1 m 

— In — ^ + = 

12 M 2 3 U 2 s 



^ + ^n^ (s ,m H ) 

(«)/ 



- 2 



f3(m 2 J S ) + 1 
s 2 J P(m%/s) ^(3(m 2 / S )-l 



In 



(68) 



5.2 The standard model case 

Quite often Eq.(55) is presented without additional theoretical support. The correct statement is that 
a is defined from Thompson scattering at zero momentum transfer. At two loops there are four classes of 
diagrams contributing to the process: 

I Irreducible two-loop vertices and wave-function factors, product of one-loop corrected vertices with 
one-loop wave-function factors; 

II one- loop vacuum polarization ® one- loop vertices or one- loop wave- function factors; 

III irreducible two-loop AA, AZ, A<p° transitions; 

IV reducible two-loop AA, AZ, Acj) transitions. 

The various contributions are depicted in Figs. 7-10, where we have made a distinction between process- 
dependent and universal corrections. 

Using QraphShot we have generated the whole set of corrections, I - IV, at p 2 = including all special 
vertices (see appendices A - C of I) .By algebraic methods, i.e. full reduction of tensor structures without 
using the explicit expressions for the scalar integrals, we have verified that 

• the non-vanishing contribution originates from III and IV only and, within these terms, only the re- 
ducible and irreducible AA transition survives. For this result the role ofT (see Sect. 8 of II) is vital. 

To present results in the full standard model we introduce auxiliary quantities, 



M 2 m 2 11 R A- 1 

i: Xi = ^ f\(M) = ---A uv (M 2 ), L^ln^, (69) 



and obtain the following results: 
— bosonic part 

r(i) 



n bo S = 6 f-i(M) + 8 x w - (3 + 4 x w ) 0{x w ) L p {x w ), (70) 
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(a) 



(b) 



(<0 




(d) 



TO 



Figure 7: Two-loop process-dependent corrections to Rutherford scattering not involving propagator corrections. 



leptonic part 



np 3 =-|/ii(M)-y^i + ^ ln:ri + ^(l + 2«O0(«OM*j), 



top - bottom contribution 



n 



(i) 



4 °/±i(M)-^+4 £ Qj[lnar / -z / + (l + 2^) J 8(« / )£ /J (* / )' 



27 



with Zi = mf/s. In the limit s — > the vacuum polarization becomes 

< o ) s (0) = 6/i 1 (M) ! 

n«(o) = -|/l 1 (M) + | + |in^ ! ng)( ) = -^/i 1 (M) + p + 4 2 



3 '93 " « " ' 9 

First we consider fcrmion mass rcnormalization, obtaining 



f=b,t 



W / = ^-( 1 + 2 l£ 



2 ^4 

2 £ 



with renormalization constants given in Sect. 5 of II. Consider the fermionic part of IlW relative 
fcrmion generation (yi,l,t and b) and perform fermion mass renormalization; we obtain 



n 



(i) 

fer 



n (1) +^_An (1) 

11 fcrm + o c AA11 fcrn 

7T e 



When we add the two-loop result we obtain 



_9_tt(D 



16 7T 



2^ 



(16tt 2 ) 5 



n (2) = one loop + 9-r \R {2) e~ 2 + e" 1 + H fin 



(71) 
(72) 

(73) 
(74) 

(75) 
to one 

(76) 
(77) 
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(1.) 



(i) 



Figure 8: Mixed corrections to two-loop fine-structure-constant renormalization equation involving propagator, vertex 
and wave-function corrections. Solid lines stand for a photon, a Z boson or a Higgs-Kibble scalar. 
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Figure 9: Process-independent irreducible corrections contributing to the two-loop renormalization equation for the 
fine-structure constant. 



The two residues are given by 



R 



(i) 



11 

256 
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~ 128 1 
+ 64^ 



A,. 



407 9 
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128 



11 

256' 



Cq x w 
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6912' 
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x H x H 



9 / x 2 1 
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128 \ x H x„ 



512 



4c fl 



ik XB + T^ XT+rlMLf3M 



18c^ Xw -\- 8 9x s J 128 



11 3 

768 + 64 



c e 4 — (1 +®%w) 

Xh 



Xh Xh X h 



9x H — 13 + 36x w x L 
36 „ 

\-9x H +4 

X M 



U4x w -± ) + —x 



9 

4— 

■ ' ' H 



x 2 
4_L 



(78) 



(79) 



Therefore mass renormalization has removed all logarithms in the residue of the simple ultraviolet pole for 
the fermionic part while a non-local residue remains in the bosonic part. 

If we work in the 't Hooft - Feynman gauge a simple procedure of W mass renormalization is not enough 
to get rid of logarithmic residues in the bosonic component and the reason is that in a bosonic loop we may 
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have three different fields, the W, the <p an d the charged ghosts and only one mass is available, M. 

The situation is illustrated in Fig. 11 where the cross denotes insertion of a counterterm 5Z M ; the latter 
is fixed to remove the ultraviolet pole in the W self-energy and one easily verifies that the total in the second 
and third line of Fig. 11 (0 and X self-energies, respectively) is not ultraviolet finite. 

The procedure has to be changed if we want to make the result in the bosonic sector as similar as possible 
to the one in the fermionic sector. Also with this goal in mind we have introduced the gauge in sect. 3 
of II. 

Collecting all diagrams, renormalizing the W mass and inserting the solution for the renormalization 
constants in the gauge we find a convenient expression for the bosonic, one-loop, AA self-energy: 



n bo S - 6 f-i( M ) + 6 + 8x w - p- x {x w )L {x w ){2> - 8x w - 16x 2 w ) + 
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with a correction given by 
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h -x H + -x L + -(x T + x B ) + 3 
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(82) 



Including both components and taking into account the additional contribution arising from renormalization 
we finally get residues for the ultraviolet poles which show the expected properties: 
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Figure 10: Process-independent reducible corrections contributing to the two-loop renormalization equation for the 
fine-structure constant. 
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Eq.(83) shows complete cancellation of poles with a logarithmic residue; furthermore the two residues in 
Eq.(83) are scale independent and cancel in the difference H(p 2 ) — 11(0). 

A final comment concerns the Z-photon transition which is not zero, at p 2 = 0, in any gauge where £ 7^ 1 
even after the Ti rc-diagonalization procedure. However, in our case, the non-zero result shows up only due 
to a different renormalization of the two bare gauge parameters and it is, therefore, of O (g 4 ) ; it can be 
absorbed into T2 which does not modify our result for II since there are no Independent terms in the AA 
transition (only appears). 

It is important to stress that in computing II one-loop field and charge counterterms are irrelevant and 
only mass and gauge parameter renormalization matters. 

Eq.(66) is not yet a true renormalization equation since the latter should contain the physical electron 
mass m e and not the intermediate parameter m R but the relation between the two is ultraviolet finite. All 
of this is telling us that a renormalization equation has the structure 

Pphys = / Q , Pbarc^ , (84) 

where the residue of the ultraviolet poles must be local. A prediction, 

O Q , Pbaro^ = O(Pphys), (85) 

gives a finite quantity that can be computed in terms of some input parameter set. 

Note that we could follow an altogether different approach, avoiding the explicit introduction of coun- 
terterms; instead of the chain bare — > renormalized — > physical we simply drop the second step, bare — ► 
physical. 



6 Running of a beyond one-loop: issues in perspective 

In the spirit of effective, running, couplings the role played by the running of a has been crucial in the 
development of precision tests of the SM. This is closely related to an hidden thought, universal corrections are 
the important ingredient, non-universal ones should be made as small as possible and represent a (somehow 
unnecessary) complication to be left for a true expert. Then, universal corrections must be linked to a 
set of pseudo-observables and all data should be presented in a way that resembles, as closely as possible, 
the language of pseudo-observables. The crucial point is that this language is intrinsically related to the 
procedure of resummation, but the latter comes into conflict with gauge invariance. The whole scheme 
received further boost from precision physics around the Z mass scale where it is relatively easy to perform 
a discrimination between relevant and (almost) irrelevant terms (boxes, for instance, are of little use) paying 
a very little price to gauge invariance. 

But what about a definition of the running of a at an arbitrary scale? What about large energies, where 
Sudakov logarithms [37] come into play? Many words have been spent in order to describe the problem and 
to derive a reasonable solution. The problem is far from trivial since large fcrmionic logarithms must be 
resummed, according to renormalization group arguments. 

One idea is to import from QCD the concept of M S couplings [38]. QCD is a theory without any obvious 
subtraction point and one defines an MS coupling where, working in dimensional regularization, poles are 
thrown away and the arbitrary unit of mass is promoted to become the relevant scale of the problem. Then, 
in QED, we take advantage of the fact that at p 2 = the vacuum polarization is gauge invariant, something 
that can be proved, to all orders, by using Nielsen identities [21]. Therefore, the idea is to express theoretical 
predictions by means of a resummed, MS, coupling [39]. This choice, is not immune from criticism: following 
common wisdom, it is more physical to use an effective charge as determined from experiment to define the 
fundamental coupling. Nevertheless, we introduce 
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Figure 11: W mass counterterm insertion for charged transitions in the 't Hooft - Feynman gauge. While the WW 
one is ultraviolet finite the same is not true for <j><j) and ghost-ghost transitions. 

The gauge parameter independence, at the basis of the MS definition, deserves an additional comment. So 
long as we are considering the one loop case in the i?j gauge, we have 
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In any gauge where £3,3(0) ^ (for a definition of the LQ basis see Sect. 9 of I) one has to take into account 
that the correct factor is 



n^(o)-^4>) = -5 r " 



- - A UV (M 2 



(88) 



gauge invariant by inspection. In our case we have 



1W0, {£}) - Wo, {£}) - n QQ (o, {£ = 1}) - S3q( o, {e = 1}) = n QQ (o, {£ = i», (89) 

due to our T prescription. 

The definition of H MS requires some explanation; we introduce the ultraviolet decomposition 



n QQ (o) = n«(o ; fc)^(Af 2 )+ ]T nf Q (o ; k) f|(m 2 ), 



(90) 
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Figure 12: The recipe for renormalizing mass dependent ultraviolet poles in the charged sector. 
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and A„ v = 7 + ln7r + lnM 2 /^ 2 . Furthermore, for 11^(0) we use 
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where the quantity within brackets is related to Aajj ad (s) and the light quark component, ReII lq (s), is 
computed in perturbation theory. The MS prescription amounts to project n QQ (0) according to 



MS 



0, V- wn e = Q, TWA, 



In- 



M 2 



(93) 



'MS 



22 



Conventionally we shall use s-^-g — s. Our results are shown in Tab. 2. In our results one-loop always 
includes the non-perturbative hadronic part; therefore the difference between one-loop and two-loop is fully 
perturbative. We observe mild variations in the range of energies from M z to 500 GeV induced by a ±5 GeV 
change in m t an by varying M H between 150 GeV and 300 GeV. 

A possible solution of the puzzle of resummation could be: do the calculation in an arbitray gauge, select 
a gauge parameter independent part of self-energy, perform resummation while leaving the rest to ensure the 
same independence when combined with vertices and boxes. The obvious criticism for this procedure is that, 
although respecting gauge invariance, it violates the criterion of uniqueness since any gauge independent 
quantity can be moved back and forward from the resummcd part to the non resummed one [40]. Let us 
accordingly define 

1 + Aa(«); (94) 



a(s) 

our result are shown in Tab. 3. The hadronic, non-perturbative, part has been computed with the help of 
the routine THADR5 by F. Jegerlehner [41]; the full calculation has been performed with CoopBack [42]. 

Are there ambiguities in the definition of the MS parameter? There is the question of decoupling of heavy 
degrees of freedom, for instance Veltman theorem is violated. One option is to introduce an MS scheme 
with decoupling of high degrees of freedom, for instance following the recipe of [43] according to which the 
In m t /M z terms in ReIT 7Z (M^) are subtracted for m t > M z . More generally, the idea is to subtract all 
contributions that involve particles heavier than the Z boson and do not decouple. Unfortunately, this is 
fully equivalent to shift portions of the result between resummed and non-resummed components. 

The behavior of Rea _1 (s) for y/s = 200(500) GeV and different values of mt,M H is given in Tab. 4. In 
Tabs. 5-6 we compare the running of the e.m. coupling constant between the MS'-scheme and our option. 
Sizable differences are present for a light H boson mass. 

To summarize we may say that ctjfg(s) is the only definition of the running of the e.m. coupling constant 
- beyond one loop - which is anchored to a rigorous, formal, basis; nevertheless, this definition is far from 
physical intuition, a pseudo-observable should always reflect the image of a (indirectly-) measured quantity. 

If we disregard the possibility of importing MS quantities into the electroweak theory only two choices 
are left: to present the full calculation for an arbitrary process, without an attempt to introduce universal 
corrections (modifications of the structure of the photon propagator, like in the Uhcling effect) w.r.t. process- 
dependent corrections, or to introduce universal corrections according to some convention. The latter choice 
requires agreement in the community. 

Any convention should be judged by the quality of the results that can be predicted. Giving a definition of 
the running e.m. coupling constant has a link to the idea of introducing some improved Born approximation, 
a concept which is also questionable at very high energies. It is for this reason that we do not give any 
special emphasis to our numerical results (if not for the fact they are there, with some degree of novelty). 
We observe sizable corrections for a relatively low Higgs mass scenario and energies well above M z ; the 
running of a(s) is almost doubled for a Higgs mass below 300 Gev and y/s = 200 Gev or higher. At the same 
time we have studied the behavior for large values of the Higgs boson mass and y/s < M H , shown in Tab. 7. 

Our definition is also affected by the presence of any threshold, e.g. the 2m t normal threshold (sub- 
leading Landau singularity) which is enhanced at two loops by the coupling with a Higgs boson. Finally and 
always at two loops we observe a strong correlation between the range of variability of the top quark and 
the Higgs boson mass. 

A final comment should be devoted to asymptotic limits: for several years important results have been 
obtained for leading and sub-leading results in various heavy mass limits, noticeably M H or m t — > oo; 
unfortunately, infinity is not always around the corner, and leading effects are masked by - large - constant 
terms, as already observed by [38]. 

The observation that asymptopia may affect upper bounds but not central values for (pseudo-)observables, 
should not be confused with a criticism to important evolutions in our recent history; only additional progress 
will help in clarifying this issue and we consider our calculation, at non-zero momentum, as a tiny step in 
this direction. 
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7 Complex poles: the paradigm 



To write additional renormalization equations we also need experimental values for gauge boson masses. 
For the W and Z bosons the IPS is defined in terms of pseudo-observables; at first, on-shell quantities are 
derived by fitting the experimental lineshapes with 

where N is an irrelevant (for our purposes) normalization constant. Secondly, we define pseudo-observables 



M P — M os cosip, T P = r os sin?/', tp = arctan — — . 



(96) 



which are inserted in the IPS. At one loop we can use directly the on-shell masses which are related to 
the zero of the real part of the inverse propagator. Beyond one loop this would show a clash with gauge 
invariance since only the complex poles 

s v = nl-i'jvfMv (97) 

do not depend, to all orders, on gauge parameters. As a consequence, renormalization equations change 
their structure. There is also a change of perspective with respect to old one-loop calculations. There one 
considered the on-shell masses as input parameters independent of complex poles and derived the latter in 
terms of the former [44]. 

Here the situation changes, renormalization equations are written for real, renormalized, parameters and 
solved in terms of (among other things) experimental complex poles. When we construct a propagator 
from an IPS that contains its complex pole, say s v , we are left to consider a consistency relation between 
theoretical and experimental values of -y v . If instead, we derive s w from an IPS that contains s z , this is 
a prediction for the full W complex pole. Note that there is a conceptual difference with calculations that 
relate MS and pole masses of gauge bosons. 

Two type of relations are known in the literature [45] : complex poles in terms of bare (or MS) masses 
and their inverse, where the MS'-masses are expressed in terms of complex poles. 

Furthermore, consistently with with an order-by-order renormalization procedure, renormalized masses 
in loops and in vertices will be replaced with their real solutions of the renormalized equations, truncated to 
the requested order. From this point of view there is no problem in our approach with cutting-equations and 
unitarity. This scheme can be used for computing pseudo-observables, including decay widths, but additional 
complications arise when we consider processes 2 — > n with unstable particles in any annihilation channel. 

Alternatively, one could use Dyson resummed (dressed) propagators, 

= i^ro < 98 » 

also in loops, say two-loop resummed propagators in tree diagrams, one loop resummed in one-loop dia- 
grams, tree in two- loop diagrams [46] . Dressed propagators satisfy the Kallcn - Lehmann representation and 
processes with external unstable particle should not be considered. This recipe requires that only skeleton 
diagrams are included (no insertion of self-energy sub-loops) and, once again, cutting-equations and unitarity 
of the S -matrix can be proven; we postpone a full discussion to Section 10 and simply mention that 

A+(P 2 ) = 0( P0 ) [A v (p 2 )] 2 2zRcE vv (p 2 ), (99) 
while, for a stable particle, the pole term shows up as 

A+(p 2 ) = 0(po) \a v ( P 2 )] 2 2iRcZ vv (p 2 ) + 2in 5(p 2 + m 2 v ). (100) 
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8 Renormalization equations and their solutions 



Renormalization with complex poles has more in it than the content of Eq.(97) and is not confined to 
prescribe a fixed width for unstable particles; it allows, al least in principle, for an elegant treatment of 
radiative corrections via effective, complex, couplings. The corresponding formulation, however, cannot be 
naively extended beyond the fermion loop approximation [44]; this is due, once again, to gauge parameter 
independence. We formulate the next renormalization equation in close resemblance with the language of 
effective couplings and will perform the proper expansions at the end. At the same time we formulate 
different choices for IPS. 

To proceed further, we also need residual functions defined according to 

Z b (s) = Z 3q (s)+F b (s), B = W 7 Z, andtf, (101) 

and discuss solutions of the renormalization equations for different IPS. One of the ingredients in our equa- 
tions is given by dressed propagators; for the W boson it has been defined in Eq.(95) of I, whereas propagators 
and transitions in the neutral sectors are given in Eqs.(106)-(108) of I. 

Furthermore, in Sect. 6 of I we have defined vector boson transitions to all orders, e.g. 

oo , 2 \ n 

D AA = s 2 e U QQ . ext p 2 = s 2 J2 ( ) 4« ; «t ( 102 ) 

n=l ^ ' 

etc, distinguishing the 6 dependence originating from external legs and the one introduced by internal legs. 
In this section we will drop the suffix ext. 



8.1 Running parameters 

As a consequence of introducing higher order corrections the coupling constant g will evolve with the 
scale according to 

^^-T^n'l-w-^ngw. (io» 

The running of e 2 = g 2 Sg is controlled by 

r ol 

e 2 (s) l-—U R (s) =4ira (104) 

L 4 7T 

while the running of the weak-mixing angle is defined according to 

s%,) - SM. (105) 

Eqs.(103)-(105) still contain renormalized parameters and in the following sections we will show how to 
replace renormalized quantities with pseudo-observables of some IPS. 

The fact that we use the relation M — M/cg, valid for bare and renormalized parameters, should not 
bring confusion; as shown in the following section the renormalization equation for eg depends on the IPS. 
In Section 8 we will introduce additional, secondary, running parameters. 

Obviously our running parameters and the MS ones are different objects and only the former have a 
physical interpretation, while the latter are nothing more than a convenient way of expressing the bare 
parameters of a renormalizable theory. 



8.2 General structure of self-energies 

In this section we clarify the issue of gauge parameter independence, order by order, of self-energies 
evaluated at their complex pole. We now examine more carefully a two-point function to all orders in 
perturbation theory, 

oo 

EvvM = E S ^M<? 2 "- (106) 

n=2 
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All one-loop self-energies corresponding to physical particles are gauge-parameter independent when put on 
their, bare or renormalized, mass-shell and coincide with the corresponding £ = 1 expression, i.e. 



^From arguments based on Nielsen identities we know that 

d 



— £ vv (s P ,£) = 0, 



s P — My + T, vv (s P ) = 0. 
To proceed further, we write a decomposition into independent and ^-dependent parts, 

and use the relation between bare (gauge parameter independent) mass and complex pole, 

Ml = s P +g 2 YSl + 3 4 [E« ;J ( Sp )S« ;€ ( Sp ,0 - - 4 2 ). c ( Sp ,0] +0 (ff 6 ) 

to derive, as a consequence of Eq.(108) and of the fact that M is a bare quantity, 

4v;€(«P.O = S|rviV(*p)*vv(*P,0, 

etc. As a result we can prove that all ^-dependent parts cancel, 



'S vv (s P ) — S vv . j(s P ) g 2 ™. 
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(107) 



(108) 



(109) 
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However, this example shows how an all-order relation should be carefully interpreted while working at some 
fixed order. 



8.3 Outline of the calculation 

Here we introduce renormalization equations and their solutions. One of our renormalization equations 
will always be of the type 

2 

s v= M v- t|-^S vv ( Sv ,M 2 ). (113) 

ID 7T 

We have three options in using Eq.(113) which, we repeat, is gauge parameter independent if the self-energy 
is considered to all orders. 

I) working at O (g 4 ) we use Eq.(113) as it stands, i.e. in £ vv we keep p 2 = —s v ; the result of Section 8.2 

guarantees absence of gauge violating terms of O (<? 4 ) (violation is due to the missing S^.^(s P ,£) term); 

II) we replace T, vv with Yi vvu which is the correct recipe but requires working in a gauge with arbitrary 

(renormalized) gauge parameters; however, up to O (g 4 ) we can use an explicit calculation for Y^l- I {s) and 
for 3>vv(s,£) (e.g. from Ref. [10]). Using Eq.(lll) we immediately derive 



; l( s v) — ^vv{ s V > £ — 1) ^vv ; l{ s v) $vv(Si/ , £ — 1). 

Ill) we expand, as done in ref. [45] 
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where the suffix p denotes derivation, 



£&(M*,M*) - ^(M'^MD^.piM^M 2 ) (115) 
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and invert, obtaining M 2 in terms of fj, 2 — Res v . The combination within brackets in Eq.(115) is gauge 
invariant (as shown by explicit calculations [45]) while Y,vl(M 2 , M 2 ) is not. 

Whenever M v can be reconstructed from other pseudo-observables of the IPS not involving s™ P (but in- 
volving other experimental complex poles) , 
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with coefficients 
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(119) 
(120) 



Note, however, that for expanding a function /(s exp ) around Res exp one has to assume Ims cxp to be O (g 2 ), 
where g 2 is expressed in terms of pseudo-observables of the same IPS. This is, for instance, needed in deriving 
M from an IPS containing Sz Xp . 

8.4 Notations 

Residual functions F B (s) for B = W,Z and H are defined in Eq.(lOl). All functions are expanded up to 
second order, e.g. 



16 7T 



2 Fw 



(2) 



Furthermore, we introduce 



Si n) (s) = # (0) - Re £ (s) + Rc £<J (a) , 



(0), Fi" ) ( S )=ReE^(s)-Fi" ) (0), 

4")( S ) = S^( S )-4" )( S ), 

(n) (o), 



(121) 

(122) 

(123) 
(124) 

(125) 



where s z = [i 2 z — i 7 Z /i z and S H = n H — i ^ H ji H are the (input) Z and H boson complex poles. The LQ 
decomposition has been introduced in Sect. 6 of I (Eqs.(123)-(125)) and, here, we drop the suffix ext. 

8.5 Input parameter set I: a, G F and M w 

The starting point for our analysis, as for the solution of rcnormalization equations, is to fix an IPS. 
Our first choice is the {a, G F , M w } IPS; although wc still call it {a, G Fl M w }, it is clear from the previous 
discussion that something slightly different is meant. Wc use a, G F and [i w and predict, among other things, 
7 W which, in turn, can be compared with the measured OS T w . We begin with two equations 
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M 2 — 9 



16 7T 
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S 3Q (s w ) + F w (s w ) 



(126) 
(127) 



The (finite) mass counterterm of Eq.(127) is to be contrasted with the conventional mass renormalization 
where ReT, ww (M 2 ) is used. 
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We look for a solution with the following form: 

G 



f = SG/i 2 w 



1 + E c » 



n=l 



M 2 = M 2 



C \ ™ 



(128) 



where G is the process independent coupling constant of Eq.(40). A straightforward calculation shows that 



C 9 {1) = \f£\s w ) C g {2) = C 2 (l) + -^ 2 w F^{s w ) C M (l) = ^Re^> w (sw) 
C M (2) = C 2 M (1) + i Re - F«(0) £&l(^)' 
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Note that there is a special combination of renormalized parameters, M 2 /g 2 , which enters into the W 
propagator; using Eq.(129) we obtain 



l + ^J#>(0) + ^ji a >(C>) 



(130) 



m 2 
IF 

We start with the running of g, 
- Renormalized running of g: 

for this IPS the complete renormalization of the coupling constant g is obtained after inserting Eq.(129) into 
Eq.(103), 



g 2 (s) SGfil 16^ 2 M 2 W 
To proceed further, we define the 
- Renormalized running of s 2 ,: 
the renormalization equation for s 2 , is 
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6g™=i&n™{s) + F™(s w ). (131) 
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Using the first of Eq.(128) we obtain a solution given by 
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( 5 s («) = Fi")( Sw ) + 25 2 M 2 w 4" Q ) (0), 



Note that we have a residual dependence on s 2 , in (5s( 2 ); here s$ must be set to the lowest order value s 2 . 
- Renormalized W propagator: 



(133) 



for the W propagator we factorize a g 2 , insert the solution and write its inverse as 



g 2 A w (s) 
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,g 2 (s) 8G 16 7T 2 
Using Eq.(131) the same expression can be rewritten as 
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where the remainders are: 



Rw\sw) = ImS^(s w ) - Hwlw II^(s w ) 



-(") 



(«), 



(136) 



The complex zero of this expression is the theoretical prediction for the complex pole of the W boson. The 
real part will differ from ^ only in higher orders, the difference being proportional to 



(!) f.fi * 
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g 2 Re e( v | v (s w ) — E^^(/x 
which is (g 2 7 W ); the solution for the imaginary part is 
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(138) 



with coefficients 

71 = ImS^(/i^), 

72 =ImE^(^) [Re^'^J-^'tO)] [lmi4 2) (/4) - Imi^V^) ReE^ iP (^)] , (139) 

where the suffix p denotes derivation. We have one consistency condition obtained by comparing the derived 
width of Eq.(139) with the experimental input j w . The goodness of the comparison is a precision test of 
the standard model. 

Furthermore, the parameter controlling perturbative (non-resummed) expansion is G F \£ w and we derive, 



G=^fl-^^^+[2(4 1) ) 2 -44 1) C 3 (l)-4 2) l (^%) }■ 
V2 >- 2V2tt 2 L /j, 2 w 9W J V2V2tt 2 / J 

In other words, we can go from the G option of Subsection 3.4 to the G F option by replacing 
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Fw = F {2) (0) + M 2 W 6%> + 5£> fil 6%> + Re F$> (s w ) + Re (s w ) - 2 F% , (141) 
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(1), 
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and G — > G F j\[2 in all the results of this section. 

All functions appearing in the results depend also on internal masses, M etc. Therefore we always use, 
for and arbitrary / 



/W^M 2 ,...)^ 1 ^;^ 



(142) 



A last subtlety in Eq.(134) is represented by the residual s 2 , dependence of the W self-energy and of S G ; we 
use 



s 2 - s 2 



2tt 2 



(143) 



8.6 Input parameter set II: a, G F and ikf z 

Within this IPS we use the following three renormalization equations: 



G 



M 



9 



16 7T- 



F w {0) 



Air a 



1 



9 



16 7T 



n QQ (o) 



„2 2 

.9 s e 

16 7T 2 



Rcn QQ (s z ) - Mz c 2 + M : 



16 7T ; 



• ReE zz (s z ) 



16 7T 2 



2 2 

= .9 s 



Re S ' z(Sz) =0. 



(144) 
(145) 
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Within this IPS we look for a formal solution of the renormalization equations which improves upon fixed 
order perturbative expansion. Once this solution is obtained we will discuss the necessary steps to reinstall 
gauge parameter independence. The improved lowest order solution for s 2 , is defined by 



2 Giil 



4 n a(s) ~ e 2 (s) , a z = Rea(s z ). 



A solution of Eqs.(144)-(145) is written according to the following expansion: 



g 2 s 2 = 4ira z 



71=1 



M 2 =y 2 z c 2 



l + Y,C M {n)a n {s z ) 



n=l 



2 -2 

si = -s- 



1 + C s (n)a n (s z ) 



n=l 



(146) 



(147) 



(148) 



To explain once again our procedure, we may say the following: despite its intrinsic simplicity, s 2 of Eq.(146) 
has problems with gauge parameter independence and we can a) expand a z up to O (a 2 ) (improving, in any 
case, upon one-loop results), b) express a z in terms of 0-^(0) and c) introduce a (non- unique) version of a z 
where the resummation is performed in terms of a gauge parameter independent choice of II H with the rest 
expanded up to second oder. To study gauge boson complex poles it is important to have a resummation of 
large fermion logarithms and the use of II H ; QED suffices. 

A solution to Eqs.(144)-(145) is provided in terms of the function S^"' of Eq.(122) and of 



The results are given by 

^( 1 ) = i ^Re[^nW( Sz ) ^ 



(149) 



L M1 
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(c 2 — S 2 ) \1 

2n^ ) ( Sz ) + ni 1 2( Sz ) 



+ 
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C M (l) = i ^Re{-^[3 Sz n«( 



Sz) + ^2 ^33 ( s z) 



(c 2 -s 2 )/4 



(150) 



Furthermore, we introduce 



-c 2 Si 2) + (1 - 3 s 2 + 2 s 4 ) Re sgj (s z ), 



c 2 + Re {s z [s 2 (1 - 3 s 2 + 2 s 4 ) U^>(s z ) + (1 - 4 s 2 + 5 s 4 - 2 J b ) 11$ (s z ) }, 
/4 AM< 2 > = -s 2 £<?> + (c 2 - s 2 ) Re [e$( Sz ) - s 2 E<g( Sz )~ 
Ae (n) = AgW + A*< n >, 



to obtain 

Cg(2) = 
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16 7T 2 s 4 c 2 (c 2 -s 2 ) 



, C a (2) = 
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(151) 
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(153) 
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C' s (2) = AsW + § 2 c 2 Rc { 16 s 4 [A^ 1 



2 /; 2 
+ 16s 4 ^ 

s z 



Ae«l + 16s 2 (l-3s 2 )Ai< 1 )Ae« 



+ n«( Sz ) f 8 S 2 A5« 



}• 



(154) 



C^(2) = AM< 2 ) + s 2 c 2 Rc {l6 s 4 [AM' 1 '] 2 - 16 s 4 [a#« + AM* 1 )] 2 + 16 s 2 Ag« AM* 1 ' 



+ 16^4 (l _Mz 



Ae« 



-n£>(„) [s^a.W-^^ 



} 



The solution that we have obtained will be interpreted in terms of running quantities: 

- Renormalized running of g: 
The running of g, given in Eq.(103), is fixed in IPS II by the following equation: 



1 s 2 
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At the same time the ratio M 2 / g 2 becomes 



M 2 _ s 2 c 2 nl 
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AlTCtz 16 7T 2 
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(155) 



(156) 



(157) 



- Renormalized running W mass: 

At this point we introduce a new running parameter, M(s), through the relation 

M 2 (s) s 2 c 2 n 2 z _L -(i) a z ^2), . , . 

g z (s) 47ra z 16 7T 2 64 

with F^(s) given in Eq.(123). At s = we have g 2 (0)/M 2 (0) = 8 G and, therefore we can describe the 
running of G by introducing 

- Renormalized running of G: 

G(S) 4^)' G(0) = G - (159) 

The running parameters introduced so far allow for a simple representation of the W propagator, one of the 
primary quantities in discussing renormalization of the standard model. 
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Renormalizcd W propagator: 



Using Eqs.(156)-(157) we obtain the following expression for the inverse W propagator (a factor g 2 is 
factorized) ; 



5 2 A w (s) 



where the residual term is 
1 



_s M 2 {s) 

g 2 (s) g 2 (s) 

S S w 



+ R,c R w — 
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(160) 
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(161) 



Eq.(160) shows the intrinsic beauty of the language; the W propagator is entirely written in terms of running 
couplings. Another useful combination of renormalized parameters is 
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(162) 



To go from the G option of Subsection 3.4 to the G F option we perform the following replacements: s 2 is 
now the solution of 



G F — 



ira 2 



2/4 s 2 c 2 ' 



(163) 



and everywhere we perform the following replacements: 

i4 1} (o) - 



^=^(0) 



■■^( Sz )-2s 2 c 2 ^(s z )-(2c 2 -s 2 )F^]}. 



■Re 



(164) 



Next, we consider the Z propagator; however, it is more convenient to introduce a new parameter which is 
related to the custodial SU(2)v symmetry: 

Running p-parameter: 

To write the Z propagator we introduce a Vcltman running parameter, p(s), defined by 

1 . a z 



1 + — 

p(s) ir 



4 7T S 
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2p 2 z 5p^(s) 



1 1 
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O 1A Z 



(165) 



where we have used Eq.(124). This new parameter has its own relevance insofar it allows us to simplify the 
Z propagator: 
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- Renormalizcd Z propagator: 
With these ingredients it can be shown that 

4a- i ( S ) = 



c 2 (s) M^s) 

g — 

g 2 (s) g 2 (s)p(s] 



+ R z (s), 



where c 2 (s) = 1 — e 2 (s)/g 2 (s) and where the residual term is given by 
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(167) 



the second order correction is 
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(168) 



One final comment concerns the gauge parameter independence of (finite) renormalization; Eqs.(127) and 
(145) are the definition of W and Z boson complex poles (actually their real part) and therefore £- 
independent. The second Eq.(144) involves 11(0) which is gauge parameter independent. As far as Eq.(126) is 
concerned we have to understand it as properly expanded to the requested order since only M 2 S G + E lvlv (0) 
is ^-independent. 

8.7 Including the Higgs boson 

To describe the Higgs boson lineshape we define the following quantities: 



„th 



M 2 H -iM„T f 



(169) 



To study properties of the Higgs boson a new equation is added to Eqs.(144)-(145). 

„2 



Ph = — 
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Re 



Z 3Q (S H ) + F H (S H ) 



(170) 



where S H is the experimental H boson complex pole. Expanding M 2 



n=l 



(171) 
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we obtain the following solution: 

M 2 H ^(l) = iRcE«(^) 

A C h M {2) = 1 Re \y£1{S h ) + \ y£1{S h ) \y£1{S h ) - F^(0)] }. 
The ratio /g 2 becomes 
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(172) 
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(174) 



- Rcnormalizcd running H mass: 
After introducing the running Higgs boson mass, M H (s), through 

Ml(s) _S^l 1 _m 



g 2 (s) Ana z 16-7T 



FV(s)- 



64 7T 3 s 2 



i4 2) (s), 



(175) 



with F„(s) defined in Eq.(125); we use Eqs.(156)-(173) to obtain the following expression for the inverse 
H propagator (a factor g 2 is factorized), 
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with a residual term 
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Fucrthermore, if we use input parameter set I we have a mass renormalization equation expanded as 

m 2 =, 2 [i+£ (^y n - 
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with coefficients 
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- Rcnormalizcd H propagator: 
The inverse Higgs propagator becomes 

.9 2 A H ( S )- 
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With two-loop accuracy the solution for 
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The coefficients are 
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(184) 



In Eq.(184) we have introduced the notation / p (s) = df/ds and two-loop contributions include the finite 
renormalization arising from one-loop terms, i.e. 



F« (s ; , M 2 ,...)= F« (s ; M 2 W , f H , . . .) + 
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with coefficients 



cif^, — Re E^^^jSw), 



C« = -p? w f fiW(O) - ReE^( Sw ) + E^(0), 



(186) 



etc. To give an idea of the impact of radiative corrections in Eq.(182) we present few results in Tab. 8 where 
it is evident that Re (s H — S H ) is a tiny effect while anomalous values for the imaginary part are impossible 
to accommodate within the standard model. 

8.8 A simple numerical example 

In this subsection we illustrate renormalization equations with a simple numerical example. First we 
write 

^=.(1 + ^ + ^) x = ^, (187) 
where g is the renormalizcd coupling constant and the coefficients are 



ai 



The perturbative solution is 



°G 



02 



= 4- s ^(o). 



x = 



16 ir< 



= X + X 2 (b 1 +b 2 X), X 



(188) 



(189) 



and the numerical impact of two-loop terms is illustrated in Fig. 9 where we have shown the percentage 
one-loop/Born corrections (b\ X) and two- loop/one- loop {b^X/bi) for different values of the Higgs boson 
mass; note that the large two-loop effect around M H = 200 GeV is due to an accidental cancellation which 
occurs at one loop level and that for high values of M H the ratio two- loop/one- loop clearly starts to indicate 
a questionable regime for the perturbative expansion. 
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9 Loop diagrams with dressed propagators 

In this section we describe how to reorganize perturbation theory when using dressed propagators, the 
so-called skeleton expansion: consider a simple model [46] with an interaction Lagrangian 

L=Z$(x)tf?(x). (190) 

The mass M of the $ -field and to of the -field be such that the $ -field be unstable. Let Aj be the lowest 
order propagators and Aj the one-loop dressed propagators, i.e. 

A _ A . 

(191) 



1 — A$ 1 — A^ 

etc. In fixed order perturbation theory, the <f> self-energy is given [47] in Fig. 13. Note that the imaginary 




\ 



— 

\ i 

a) skeleton b) £ insertion c) skeleton 



-- + A I---- + 

\ / \ / 

\ / N / 



Figure 13: The <f> self-energy with skeleton expansion, diagrams a) and c), and insertion of a sub-loop E$$, diagram 
b). 

part of £<^ is non-zero only for — p 2 > 9 to 2 (the three-particle cut of diagram b) in Fig. 13), if to <C M. 
When we use dressed propagators only diagrams a) and c) are retained in Fig. 13 (for two-loop accuracy) 
but in a) we use A$ with one-loop accuracy: 
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where we assume p 2 < 0. Since the complex 4> pole is defined by 
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we write the inverse (dressed) propagator as 
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(194) 



expand in g as if we were in a gauge theory with problems of gauge parameter independence and obtain 
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(195) 



where Q 2 — 1 — 4to 2 /s m . Note that there is an interplay between using dressed propagators for all internal 
lines of a diagram and combinatorial factors for diagrams with and without dressed propagators. Note also 
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that the poles in the q° complex plane remain in the same quadrants as in the Feynman prescription and 
Wick rotation can be carried out, as usual. Evaluation of diagrams with complex masses does not pose 
a serious problem; in the analytical approach one should, however, pay the due attention to splitting of 
logarithms. Consider a Bq function, 



5 b 2 ; M 1 ,M 2 ) = A l 



X(x) = -p 2 x 2 + (p 2 + Ml - M 2 ) x + M 2 -iS, 



(5 — > 0+ ) where one usually writes 



In — — — ln( — i 5) + hi(x — x—) + \n(x — %+). 



Since lm.x(x) does not change sign with x € [0, 1] the correct recipe for M 2 — m 2 — i mj is 



In 



X(x) 



In \p 2 | + \n(x - x-) + 0(\p 2 ) lnA(a; + - x) + r)(-X- , Xx + ) 



(196) 



(197) 



(198) 



Where rj is the 't Hooft - Veltman function. In the numerical treatment, instead, no splitting is performed 
and no special care is needed. This is specially true for higher legs one-loop functions and for two-loop 
functions. 

A t -channel propagator deserves some additional comment: one should not confuse the position of the 
pole which is always at /i 2 — i fi 7 with the fact that a dressed propagator function is real in the t -channel. 
Therefore, using one-loop diagrams with one-loop dressed $ propagators is equivalent, to O (g 4 ) , to introduce 



A$(s M ) 



A (m 2 ) 




+ 



+ 



A^(m 2 ) 



A (m 2 ) 



pole 




A$(s M ) 



\ / 



A (m 2 ) 



Figure 14: Diagram b) of Fig. 13 with one-loop dressed $ propagators is equivalent, up to O (<? 4 ) , to the sum of 
three diagrams with lowest order propagators with the $ mass replaced with the $ complex pole. The Z po i c vertex 
is given in Eq.(199) 

the sum of the three diagrams of Fig. 14 where $ propagators are at lowest order but with complex mass 
s M and where the vertex Z po \ c is defined by 

9 2 

Z vo \ c = ttt^ B o (s M ; m,m). (199) 

16 1T Z 

In the following sections we will discuss an extension of this simple scalar theory. 



10 Unitarity, gauge parameter independence and WST identities 

The simple scalar model of Section 9 is not adequate for describing the complexity of a gauge theory. A 
critical question is how to construct a scheme for a gauge theory with unstable particles which allows to deal 
with the calculation of physical processes at one and two loops. This scheme must satisfy a certain number 
of requisites, in particular wc seek for a scheme that 
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a) respects the unitarity of the S— matrix; 



b) gives results that are gauge-parameter independent; 

c) satisfies the whole set of WST identities. 

Resummation will be part of any scheme, a fact that introduces additional subtleties if a) — c) are to be 
respected. Consider in more details the definition of dressed propagator: we consider a skeleton expansion 
of the self-energy £ with propagators that are resummed up to O (n) and define 



A("V) = A(°V) [A (0 V) - S« (p 2 , A^V) 



where the Born propagator (tensor structures are easily included) is 

A(°V) ' 



p 2 + m 2 



If it exists, we define a dressed propagator as the formal limit 

A(p 2 ) = lim £<"V)> A(p 2 ) = A<°V) [a<°V) - S ( P 2 , A(p 2 )) 



(200) 



(201) 



(202) 



which is not equivalent to a rainbow approximation and coincides with the Schwinger - Dyson solution for 
the propagator. 




Figure 15: Schwinger - Dyson equation for the self-energy (l.h.s). In the r.h.s. the red oval is the SD vertex and the 
yellow box is the dressed propagator. 



Figure 16: Schwinger-Dyson equation for a dressed propagator (yellow box); the red oval is the SD self-energy. 

Technically speaking, one should also introduce an integral equation for the four-point functions, as they 
appear in Fig. 17. However, this is well beyond the scope of the present discussion and its structure shows 
a somewhat greater complexity. In discussing unitarity an essential tool is representend by the so-called 
cutting rules. 
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Figure 17: Schwinger-Dyson equation for a dressed vertex (l.h.s); in the r.h.s. we have SD three- and four-point 
vertices (red ovals) and dressed propagators (yellow boxes). 



- Cutting rules 

We assume that Eq.(202) has a solution that obeys Kallen - Lchmann representation 

Re A(p 2 ) =ImS(j) 2 ) Up 2 + m 2 - Re £(p 2 )) 2 + (lm£(p 2 )) 
A dressed propagator, being the result of an infinite number of iterations, 

POO 



= np(-p 2 ). 



A(p 2 ) 



ds 



p 2 + s 



i8' 



(203) 



(204) 



is a formal object which is difficult to handle for all practical purposes. 

Unitarity follows if we add all possible ways in which a diagram with given topology can be cut in two. 
The shaded line separates S from £t. For a stable particle the cut line, proportional to A + , contains a pole 
term 

A + = 2i7r9(p )S(p 2 +m 2 ), (205) 

whereas there is no such contribution for an unstable particle. We express Im £ in terms of cut self-energy 
diagrams and repeat the procedure ad libidum and prove that cut unstable lines are left with no contribution, 
i.e. unstable particles contribute to the unitarity of the S— matrix via their stable decay products. The 





Figure 18: Cutting equation for dressed propagator. 



consistent use of dressed propagators gives a general scheme where unitarity is satisfied which is essentially a 
statement on the imaginary parts of the diagrams. Approximated, or truncated, schemes (e.g. resummation 
of one-loop self energies, or rainbow approximation without further resummation of the vertex functions) 
usually lead to gauge dependent results. 

- WST identities 

We assume that WST identities hold at any fixed order in perturbation theory for diagrams that contain 
bare propagators and vertices; they again form dressed propagators and vertices when summed. We expect 
that an arbitrary truncation that preferentially resums specific topologies will lead to violations of WST 
identities. These violations, of course, are not present in the exact calculations. 
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Gauge parameter dependence 



A truncated approximation, e.g. simple resummation of two-point functions, necessarily leads to gauge 
dependent results. A convenient tool is to analyze the gauge invariance of the effective action where one can 
show that on-shell gauge dependence always occurs at higher order than the order of truncation. 

- Complex pole 

A property of the S— matrix is the complex pole A (p 2 = —s P ) = 0, which is gauge parameter independent 
as shown by a study of Nielsen identities, see Section 8.2. 

An approximate solution of the unitarity constraint is as follows: 



2ImT« = 
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(206) 



where, S = 1 + iT and where D(p 2 ) is the unknown form of the propagator. Making the approximation, 

2 



E 



dPS n 



M 



we derive lmD(p 2 ) = mr to t- A simple but, once again, approximate solution is 

D{p 2 ) = (p 2 +m 2 -imT tot y\ 



(207) 



(208) 



which is valid far from the mass shell and where the invariant mass at which the decay is evaluated is 
identified with m 2 . We can improve upon this solution by writing instead 



D(p 2 ) = (p 2 - Sp y\ 
which is equivalent to resum only the self-energy (up to some fixed order), and to use 

m 2 = s P + S(s P ), D{p 2 ) = -\s-s P - £(«) + S(s P )l 1 = - (p 2 - s P ) 



+ h.o., 



(209) 



(210) 



where higher order terms are neglected. Another way to see that Eq.(209) is an improvement of Eq.(208) is 
to observe that 

(211) 



■ - tot 2 

p + m + i p 

m 



— = 1 + i — (p* + Sp ) + h. . ~ p J + S[ 



A propagator with the correct analytical structure, p 2 — s P , will be represented with a thick dot. The 
approximation of Eq.(209) allows us to write the cutting equation of Fig. 19. One can see that using 
truncated propagators with complex poles (at the one-loop level of accuracy) is still respecting unitarity of 
the S— matrix within the approximation of Eq.(207) if the complex pole is computed from fermions only; 
however, this scheme violates gauge invariance since vertices are not included. There is a solution to this 
problem, namely replacing everywhere the (real) masses with the complex poles, couplings included; this is 
known in the literature as complex mass scheme [4] 

The complex mass scheme 

Since WST identities are algebraic relations satisfied separately by the real and the imaginary part one starts 
from WST identities with real masses, satisfied at any given order, replaces everywhere m 2 — ► s P without 
violating the invariance. In this scheme, using for instance IPS I(II), also s 2 , is expanded around a complex 
value, 



J 2 = 



1 7r a 



2 Gt 



or 2Gs z = 



ira. 



(212) 



In turns, this scheme violates unitarity i.e. we cannot identify the two sides of any cut diagram with T and 
respectively. 

To summarize, the analytical structure of the S— matrix is correctly reproduced when we use propagator 
factors p 2 + s P but unitarity of S requires more, a dressed propagator 
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p 2 + Sp 


p 2 + s P - S(p 2 ) + £(- s P ) 


analyticity 


unitarity 



Another drawback of the scheme is that all propagators for unstable particles will have the same functional 
form both in the time-like and in the space-like region while, for a dressed propagator the presence of a pole 
on the second Riemann sheet does not change the real character of the function if we are in a t— channel. 

In some sense the scheme becomes more appealing when we go beyond one loop. As a preliminary step 
we have verified that WST identities are satisfied with bare (i.e. non-dressed) propagators and vertices up 
to two loops; we may assume that they are verified order by order to all orders, 



w {1) ({r» = w {2) ({r» = • • • = o, 



'(2) 



(213) 



where {r} is a set of (off-shell) Green function and W — is the WST identity. 

Next we write the same set of WST identities but using a skeleton expansion with one-loop dressed 
propagators. Calling the scheme complex mass scheme is somehow misleading; to the requested order we 
replace everywhere m 2 with ,s P + E(sp) which is real by construction. If only one-loop is needed then 



s P everywhere (therefore justifying the name complex mass) and 

WW({T}) „ =0, 



is trivially true. Also, 



w (2) ({r» 



= o. 



We have two-loop diagrams with no self-energy insertions where m 2 = 



s P + £(s P ) and the factor 



E(p 2 ) - E( Sp ) 



p 2 + s { 



(214) 
(215) 

and one-loop diagrams where 

(216) 



expanded to first order with £ = £W. Furthermore, in vertices we use m 2 = s P in two-loop diagrams and 
m 2 = s P + E(sp) in one-loop diagrams. Expanding the factor of Eq.(216) generates two-loop diagrams with 
insertion of one-loop self-energies plus one-loop diagrams with one more propagator and a vertex proportional 
to E(sp); furthermore one-loop diagrams with m 2 dependent vertices get multiplied by E(s P ); it follows that 



w(i+2) / rpi , \ 

yv \\i J skclcton/cxpandcdy 

as a consequence of Eqs.(214)-(215). 



m 2 — sp+H(sp) 



0. 



(217) 
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In this scheme, one replaces Eq.(127) by 



= M 2 - 




^3q(s w ) + F w (s w ) 



(218) 



and organizes the perturbative expansion according to strategy outlined in this section, namely, skeleton 
expansion expanded around complex poles. Note that M is always a real quantity, v 

11 Beyond the complex mass scheme: outlook 

A possible improvement of the complex mass scheme which also respects the unitarity of the S— matrix 
can be constructed when we observe that unstable particles contribute to the unitarity of the S— matrix via 
their stable decay products, i.e. light fermions. Furthermore, at one loop, the so called fermion-loop scheme 
respects all WST identities. 

The birth of a complex pole is full of subtleties; the origin is the instability of the corresponding quantum 
state, reflected by the fact that the self-energy develops an imaginary part. Several statements that are 
usually made in this context are questionable and should be understood with the due caution. We have 
a few ingredients, a bare (renormalized) mass and the corresponding complex pole; we observe that Re s v 
and M 2 differ because of ImS vv ^ 0. A criterion of naturalness (which also should be taken with the due 
caution) requires a small difference, fully accountable in perturbation theory. 

Consider the toy model of Section 9; the usual statement that the mass M of the $ -field and m of the 
4> -field be such that the $ -field be unstable is based on the assumption of naturalness since both masses 
are not input data but should be derived by using some IPS, i.e. no assumption can be made a priori on 
Lagrangian parameters. Indeed one could introduce the following paradox: assume that m and M have been 
derived with some IPS and that they respect the condition m <C M; consider the <f> self-energy given by 
diagram a) of Fig. 13; look for a complex solution of 



where, for Res m large enough, ImS(s m ) 7^ 0. The paradox is that we end up with Res m m 2 and the 
notion of stable particle is lost. The resolution of the paradox lies in the observation that $ is unstable and 
must be removed from the in/out bases of the Hilbert space; in other words we cannot derive anything by 
naively cutting diagram a) of Fig. 13. This diagram, with a bare $ propagator, should not be there and only 
its dressed version should be considered; therefore, the imaginary part starts with the three-particle cut of 
diagram b) which enters whenever a virtual c/> appears in a diagram. The solution of Eq.(219) is real and 
is stable. 

Consider now a gauge theory: how is the complex pole affected by bosonic corrections? Assume that M w 
and M z (both renormalized parameters) have been derived by some IPS and that, as expected, M z < 2 M w . 
Due to gauge invariance the one-loop Z self-energy cannot develop an imaginary part, for Res z close to M 2 , 
in diagrams with W, <f> or FP ghosts so that bosonic corrections only modify Re s z and instability is a 
consequence of (light) fermion corrections. Even this statement should be slightly modified: bosonic states 
enter through their dressed propagators and when we cut the sequence stops only when cutting (at some 
higher order in perturbation theory) (light) fermion lines. Therefore, bosonic contributions change Im s z 
only at higher orders through their stable decay products. Only light fermions matter and this represents 
the basis for our proposal. 

Consider any unstable vector boson V and construct the corresponding one-loop self-energy; a rearrange- 
ment of the perturbative expansion is required. First we define the V dressed propagator according to Fig. 20, 
where only the one-loop fermion self-energy is resummed. Any triple gauge boson vertex is also computed 
in the fermion-loop approximation of Fig. 21. In constructing Sw at one loop we only include £^ v . What 
about the bosonic part? Consider the loop diagram in Fig. 22; cutting dressed propagators gives stable inter- 
mediate states with the result shown in the r.h.s of Fig. 22. Thus, the imaginary part of the diagram is one 
of the terms contributing to V — > 4 / in the fermion loop approximation. Clearly this fermion-loop improved 



m 2 - s m - 



16 TT 2 



S(-S m ) = 0, 



(219) 
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Figure 20: One-loop dressed propagator in the fermion loop approximation. 




+ permutations 



Figure 21: One-loop triple gauge boson vertex in the fermion loop approximation. 



tree diagram cannot respect gauge invariance by itself. Consider all fcrmion-loop improved tree diagrams 
contributing to V — ► 4/, as illustrated in Fig. 23; their sum is clearly gauge invariant and respects all WST 
identities. Using this total as the imaginary part of a self-energy we can reconstruct a self-energy which has 
the right properties wrt unitarity and gauge invariance. The reorganization of perturbation theory replaces 
the simple one-loop diagram in the VV self-energy with internal unstable bosons with a series of multi-loop 
self-energies where all internal boson lines represent dressed propagators as in Fig. 20. An example is given 
in Fig. 24. 

The next step in the calculation of the bosonic VV self-energy starts with the set of all tree diagrams 
contributing to V — ► 6 / written with the fermion loop improvement; after identifying this set with T(l — > 6) 
we can reconstruct another piece of the VV self-energy, TT^, using the relation of Fig. 18 and removing the 
cut on the dressed propagator. 

The same algorithm can be generalized to deal with 2^2 (or 2 — > n) processes (no unstable particles 
in the initial and final state) by constructing the fermion-loop improved T(2 — > 2(n)), by writing all the 
cut diagrams with the cut separating T from T> , using the relation of Fig. 18 and removing the cut on the 
dressed propagator. 

Computing loop integrals with dressed propagators is not the difficult part in this extension of the familiar 
scheme; they satisfy Eq.(204), so that we can compute loop integrals with masses s%, s n and fold the result 
with spectral functions, p(si). The intrinsic obstacle, for all practical purposes, is represented by the presence 
of multi-loops. 

12 Conclusions 

In this paper, we begin with the (two-loop) renormalized Lagrangian of the Standard Model and show how 
to construct (pseudo-)observables up to two loops. For this we have written a set of (finite) renormalization 
equations (RE) and shown how to solve them. These equations are based on the choice of some input 
parameter set; in this case we use the Fermi coupling constant G F , the fine structure constant a and complex 
poles for unstable vector bosons. The structure of our renormalization equations at two loops generalizes 
one-loop results that are well known in the literature [10]. 

At any order in perturbation theory quantities in the IPS - like G F - are experimental data points and, 
for the specific example, we do not make any prediction for the two-loop decay rate, no more than we do for 
the two-loop Z complex pole. These quantities are input for the procedure and allow for the determination of 
other observables. Of course, if we are really interested in the muon lifetime, then another data point instead 
of G F should be used and the lifetime would follow consistently. It is, however, an important observation 
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Figure 22: One-loop dressed propagator in the fermion loop approximation. 




+ permutations 



Figure 23: Fermion-loop improved tree diagrams for V — > 4 /. 

that the correction Ag of Eq.(6), which is known to be ultraviolet and infrared finite at one loop, only needs 
one- loop counterterms to show the same property at two loops. 

For the IPS specified by a, G F and the Z complex pole we have also considered a solution of the RE 
improved by resummation. Note that the transition from real to pseudo-observables involves certain as- 
sumptions but, within these assumptions, there is a well defined mathematical procedure. 

We have discussed the change of perspectives with respect to old one-loop calculations where one considers 
the on-shell masses, derived from a fit to pseudo-observables (line-shapes etc.), as input parameters and 
derives complex poles in terms of the former. Since only complex poles are gauge parameter independent to 
all orders it follows that on-shell masses loose their central position in any renormalization procedure which 
goes beyond one loop. 

We have also discussed how to generalize the concept of running couplings and presented numerical results 
for the value of of an effective e.m, coupling a at an arbitrary scale; this shows that not only the algorithmic 
part of our scheme is at work but also that its parallel implementation in a numerical code is complete [42] . 

It is an unfortunate consequence of gauge invariance that we cannot have a full extension of the running, 
resummed, couplings to the bosonic sector of the theory; in any case, even the concept of a fermionic sector 
is meaningless at two loops. 




Figure 24: Examples of the multi-loop bosonic self-energy replacement. 
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We have devoted the last three sections to discuss renormalization of a gauge theory in the presence of 
unstable particles, objects which should be removed from the asymptotic states of the Hilbert space (in/out 
bases). Although there is a way of formulating the general treatment of the theory, as shown by Veltman 
in the sixties in his seminal paper, one has to admit that all practical approaches have loopholes since they, 
necessarily, rely on truncation of the Schwinger-Dyson equations. 

By a suitable modification of the form of the solution of our renormalization equations we recover the 
complex mass scheme of Ref. (the less inconsistent and, so far, the only available scheme) and discuss its 
extension to two loops as well as a possible improvement aimed to satisfy the unitarity of the S-matrix. 
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A Tables 



M B os [GeV] 


150 


300 


500 


2tt 2 7tT> 
°G 


18.29% 


8.89 % 


-24.62% 



Table 1: The ultraviolet, infrared finite remainder for G F of Eq.(39) for different values of the on-shell Higgs mass. 
Here n% is the real part of the W boson complex pole, see Eq.(113). 



m t = 174.3 GcV 


M H = 150 GcV 










^ [GeV] 


M z 


120 


160 


200 


500 


one-loop 


128.104 


127.974 


127.839 


127.734 


127.305 


two-loop 


128.040 


127.967 


127.891 


127.831 


127.586 


% 


-0.05 


-0.01 


+0.04 


+0.08 


+0.22 


m t = 179.3 GcV 


M H = 150 GcV 










one-loop 


128.113 


127.982 


127.847 


127.742 


127.313 


two-loop 


128.048 


127.980 


127.911 


127.857 


127.636 


% 


-0.05 




+0.05 


+0.09 


+0.25 


m t = 174.3 GcV 


M H = 300 GeV 










Vi [GeV] 


M z 


120 


160 


200 


500 


one-loop 


128.104 


127.974 


127.839 


127.734 


127.305 


two-loop 


128.046 


127.921 


127.790 


127.689 


127.272 


% 


-0.05 


-0.04 


-0.04 


-0.04 


+0.03 



Table 2: a^fs) for different values of m t ,M H and M w = 80.380 Gev, M z = 91.1875 GeV; a s (M z ) = 0.11. The 
percentage effect of two-loop corrections is shown. 
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Aa 


value at V« = 200 GcV 


2L Re EW 


-0.00645(1) 


2L ImEW 


+0.00075(1) 


2L Rep-QCD 


+0.0005512(5) 


2L Imp-QCD 


+0.0001174(5) 


fin rcn 


-0.0000977- 0.0000998 i 


1L Re EW 


-0.06953 


lLImEW 


+0.00334 


non-pert 


-0.03266 


Full Rea(s) 


0.0078934(1) 


1L + non — pert Rea _1 (s) 


127.509 


Full Rearms) 


126.688(2) 



Table 3: Aa(s) of Eq.(94) at yfs = 200 GeV; M w = 80.380 Gev, M z = 91.1875 GeV, m t = 174.3 GeV, M H = 150 GeV 
and a B (M z ) = 0.11. p-QCD is perturbative QCD. Only the error for the perturbative calculation is reported. 



y/s = 200 GeV 






m t ,M H [Gev] 


1L + 2L 


2L/1L perturbative only 


169.3, 150 


126.774(3) 


14.64% 


174.3, 150 


126.688(2) 


16.27% 


179.3, 150 


126.598(3) 


17.97% 


169.3, 300 


127.300(3) 


4.35 % 


174.3, 300 


127.313(2) 


1.51% 


179.3, 300 


127.122(3) 


7.73 % 


s/s = 500 GeV 






m t ,M H [Gev] 


1L + 2L 


2L/1L perturbative only 


169.3, 150 


125.430(2) 


29.97% 


174.3, 150 


125.259(2) 


33.92 % 


179.3, 150 


125.070(2) 


38.29 % 


169.3, 300 


127.121(2) 


4.19% 


174.3, 300 


126.960(2) 


0.57% 


179.3, 300 


126.781(2) 


3.44 % 



Table 4: a(s) at y/s = 200(500) GeV; M w = 80.380 Gev, M z = 91.1875 GeV, a s (M z ) = 0.11 and different values 
of mt,M H . The ratio between two-loop and one-loop perturbative corrections is shown. Only the error for the 
perturbative calculation is reported. 



47 



•v/i [GeVl 




200 


300 


500 


1000 


Re oT^/Re or 1 (s ) 


1 


0.984 


0.981 


0.972 


0.957 




1 


0.999 


0.998 


0.997 


0.995 


Rea-^J/a^s) 


1.006 


0.991 


0.989 


0.982 


0.967 



Table 5: Comparing the running of alpha with the running of ctjjg for m t = 174.3 GeV and M H = 150 GeV. 





100(V5B") 


200 


300 


500 


1000 


Rea- 1 (s)/Rea- 1 (s ) 


1 


0.987 


0.988 


0.985 


0.977 




1 


0.998 


0.996 


0.994 


0.990 


Rea-^J/a^s) 


1.006 


0.996 


0.999 


0.998 


0.992 



Table 6: Comparing the running of alpha with the running of ctjjg for m t = 174.3 GeV and M H = 300 GeV. 
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Vi/M H [GeV] 


150 


300 


400 


450 


100 


128.806(2) (2.79%) 


128.831(3) (3.03%) 


128.382(8) (9.55%) 


127.893(13) (24.96%) 


200 


126.688(2) (16.27%) 


127.213(3) (5.99%) 


127.910(4) (6.85%) 


128.285(6) (12.34%) 



Table 7: a(s) for m t = 174.3 GeV increasing values of M H . The second entry gives the ratio between the perturbative 
two-loop contribution and the perturbative one-loop. 





7h 


M H 


r H 


300 


4 


299.96 


8.374 


300 


12 


299.87 


8.376 


500 


40 


500.17 


63.37 


500 


80 


500.42 


63.34 



Table 8: Values for s H derived from S H , Eq.(182). 



M° s [GeV] 


150 


200 


250 


300 


350 


h X (%) 


+3.31 


+0.13 


-2.30 


-4.84 


-7.85 


h 


+12.28 


+0.47 


-8.51 


-17.95 


-29.07 


b 2 X 


+0.25 


-1.31 


-1.38 


-2.58 


-9.26 


b £X(%) 


+2.06 


-277.31 


+ 16.16 


+14.35 


+31.85 



Table 9: Numerical solution for the renormalized SU(2) coupling constant, Eq.(189) corresponding to M w — 
80.380 Gev, M z = 91.1875 GeV, m t = 174.3 GeV. 
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